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ABSTRACT 

Spectra obtained with the Hubble Space Telescope GHRS are combined with high- 
resolution optical spectra and UV spectra from Copernicus to study the abundances 
and physical conditions in the diffuse interstellar clouds seen along the line of sight to 
the star 23 Ori. Multiple absorption components are present for each of several distinct 
types of gas, which are characterized by different relative abundance and depletion 
patterns and physical conditions. 

Strong low-velocity (SLV) absorption, due to cool, moderately dense neutral gas 
and representing about 92% of the total iV(H I), is seen for various neutral and singly 
ionized species at +20 km s _1 <; Vq £ +26 km s _1 . Most typically severely depleted 
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species are less depleted by factors of 2-4, compared to the "cold, dense cloud" pattern 
found, for example, in the main components toward £ Oph. For the two strongest SLV 
components, T ~ 100 K and the thermal pressure log(nnT) ~ 3.1 cm~ 3 K; we thus 
have tih ~ 10-15 cm~ 3 and a total thickness of 12-16 pc. The adopted average SLV 
electron density, n e = 0.15±0.05 cm~ 3 , implies a relatively large n e /nu ~ 0.01, and 
thus some ionization of hydrogen in these predominantly neutral components. 

Weaker low-velocity (WLV) absorption, probably largely due to warmer neutral 
gas, is seen primarily for various singly ionized species at km s _1 <^ v® <J +30 km 
s _1 . The depletions in the WLV gas are typically less severe by a factor of 2-3 than 
in the SLV gas, and are somewhat similar to the "warm cloud" pattern seen in lines 
of sight with low reddening, low mean density, and/or low molecular fraction. If T ~ 
3000 K for the WLV components, then we have log(nHT) ~ 4.7-4.8 cm _3 K, ran ~ 
15-20 cm -3 , n e ~ 0.2 cm -3 , re e /«H ~ 0.01, and a total thickness of 0.7-0.9 pc. 

Absorption from a number of singly and doubly ionized species, perhaps due to a 
radiative shock, is seen at —108 km s _1 <J v® ^ —83 km s _1 . While the depletions in 
these ionized components are uncertain due to unobserved ionization stages, aluminum 
(typically severely depleted) is probably depleted there by only a factor ~ 3, even at 
cloud velocities in excess of 100 km s -1 . The individual high- velocity components 
typically have T ~ 8000±2000 K, n e = ra H ~ 0.4-0.5 cur 3 , thermal pressure log(2ra e T) 
~ 3.7-4.0 cm _3 K, and thicknesses of order 0.1 pc. 

Weak absorption components from ionized (H II) gas are seen in C II, Mg II, 
and Si III at intermediate velocities (—43 km s _1 <J v® <J —4 km s _1 ). Broad, weak 
absorption from the higher ions S III, C IV, Si IV, and N V is centered at -5 km s" 1 ^ 
v® <^ +6 km s -1 . No obvious absorption is discerned from a circumstellar H II region 
around 23 Ori itself. 

The large range in n e (from 0.04 cm~ 3 to 0.95 cm" 3 ) derived independently 
from nine pairs of neutral and singly ionized species in the SLV gas suggests that 
additional processes besides simple photoionization and radiative recombination affect 
the ionization balance. Charge exchange with protons may reduce the abundances 
of S I, Mn I, and Fe I; dissociative recombination of CH + may help to enhance C I. 
The large apparent fractional ionization in the SLV and WLV gas may be due to an 
enhanced flux of X-rays in the Orion region, to mixing of neutral and ionized gas at 
the boundary of the Orion-Eridanus bubble, or perhaps (in part) to charge exchange 
between singly ionized atomic species and large molecules (in which case the true n e 
would be somewhat smaller). Comparisons of the SLV depletions and bh with those 
found for the strong "component B" (v® ~ — 14 km s -1 ) blend toward ( Oph hint at a 
possible relationship between depletion and local density for relatively cold interstellar 
clouds. Calcium appears to be more severely depleted in warm, low density gas than 
has generally been assumed. 

An appendix summarizes the most reliable oscillator strengths currently available 
for a number of the interstellar absorption lines analyzed in this work. 

Subject headings: atomic data — ISM: abundances — ISM: general — stars: individual 
(23 Ori) — ultraviolet: ISM 



-3- 



1. Introduction 

Determination of the physical properties of individual interstellar clouds from observations 
of absorption lines requires both high spectral resolution and access to the satellite UV. High 
spectral resolution is needed to distinguish the individual components contributing to the generally- 
complex absorption-line profiles along most lines of sight. High-resolution (FWHM = 0.3-1.2 
km s _1 ) surveys of interstellar Na I and Ca II absorption indicate, for example, that the median 
widths (FWHM) of individual interstellar components seen in Na I and the median separation 
of adjacent components, in both species, are all £ 1.2 km s" 1 (Welty, Hobbs, & Kulkarni 1994; 
Welty, Morton, & Hobbs 1996). Furthermore, those adjacent components may be characterized 
by very different abundances and physical characteristics. Many ground-state transitions for 
the dominant ionization states of the most abundant elements — and thus many of the crucial 
diagnostic lines for determining such quantities as the local electron density n e , the local hydrogen 
density nn, the temperature T, and the ambient radiation field — are found in the UV. The 
moderate resolution (FWHM ~ 15-25 km s _1 ) UV spectra obtained by the Copernicus and IUE 
satellites provided considerable information on total line-of-sight abundances and general physical 
conditions in the interstellar medium (e.g., Spitzer & Jenkins 1975; Cowie & Songaila 1986; de 
Boer, Jura, & Shull 1987; Jenkins 1987), but could not, in general, provide detailed information 
for individual interstellar clouds. Even detailed studies which applied profile- fitting techniques 
to high signal-to- noise ratio (S/N) Copernicus spectra of multiple transitions (e.g., Ferlet et al. 
1980; Martin & York 1982) did not discern all the components present along those lines of sight. 
The increased resolution (FWHM ~ 3.5 km s _1 ) and S/N (up to ~ 1000) achievable in spectra 
obtained with the Hubble Space Telescope Goddard High Resolution Spectrograph (GHRS) allow 
individual interstellar clouds to be studied in more detail, where the clouds are warm and/or 
isolated in velocity [e.g., the halo clouds studied by Spitzer & Fitzpatrick (1993, 1995) and by 
Fitzpatrick & Spitzer (1997)]. For the colder clouds and more complex lines of sight commonly 
found in the Galactic disk, however, even the GHRS is not able to discern and characterize most 
of the detailed component structure present. For example, the two component groups studied by 
Savage, Cardelli, & Sofia (1992) in GHRS spectra of ( Oph are found to be composed of at least 
4 and 8 components, respectively, in high-resolution Na I and Ca II spectra (Welty et al. 1994, 
1996; Barlow et al. 1995). If this underlying complexity is not incorporated into analyses of the 
UV spectra, the derived abundances and physical properties must be viewed as averages over 
several potentially quite disparate parcels of gas, and column densities obtained from moderately 
saturated lines may be systematically somewhat underestimated. 

The Orion region has provided a fruitful laboratory for the investigation of various interstellar 
processes, even apart from studies of the background giant molecular cloud. High-resolution 
surveys of optical absorption lines of Na I and Ca II have revealed very complex interstellar 
component structure in the neutral gas at relatively low LSR velocities toward many of the stars 
in Orion (Hobbs 1969; Marschall & Hobbs 1973; Welty et al. 1994, 1996). UV spectra of a number 
of stars in the region have revealed in addition absorption from ionized gas at intermediate (—80 
km s _1 wlsr ^ — 20 km s _1 ) and high velocities (—120 km s^ 1 ^ «lsr ^ —80 km s _1 ) (Cohn & 
York 1977; Cowie, Songaila, & York 1979; Trapero et al. 1996). While some of the components 
arise in relatively nearby gas within about 70 pc (Frisch, Sembach, & York 1990), much of the 
complex structure is likely to be due to interactions between the young, luminous Orion stars (and 
earlier supernovae) and the ambient ISM. Maps of the H I 21 cm emission in the Orion-Eridanus 
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region exhibit numerous filamentary structures, whose systematic behavior with velocity strongly 
suggests that some of the filaments trace the boundary of an expanding shell (Heiles 1976; Brown, 
Hartmann, & Burton 1995). Similar filamentary structures seen in Ha and [N II] emission suggest 
that the inner edge of the shell is ionized, presumably by UV radiation from the Orion OBI 
association stars (Sivan 1974; Reynolds & Ogden 1979). Enhanced soft X-ray emission at 1 /± 
and/or keV, in the region bounded by the Ha filaments, suggests that the shell encloses a 
superbubble of hot gas [log(T) ~ 6.2 K] produced by supernovae and stellar winds over the past 
several million years (Burrows et al. 1993; Snowden et al. 1995). Comparisons of the 21 cm data, 
the soft X-ray maps, the IRAS 100 /um maps, and optical interstellar absorption-line data toward 
stars in the region have yielded information on the relative locations and distances of some of the 
neutral and ionized material (Guo et al. 1995). 

This study of the line of sight to 23 Ori is the first in a series which will combine high- 
resolution optical spectra with GHRS spectra to discern and investigate the properties of individual 
interstellar clouds in diverse Galactic environments. The high-resolution optical spectra reveal 
much of the underlying component structure that is unresolved in the UV spectra, enabling the 
latter to be interpreted more accurately. In Section ||, we describe the new high-resolution optical 
and UV spectra obtained toward 23 Ori and introduce our analysis of the observed absorption-line 
profiles. In Sections || and ||, we describe in detail the profile analysis and discuss the abundances 
and physical properties derived for several "natural" groups of neutral and ionized components 
along the line of sight. In Section ||, we discuss those component groups in the context of prior 
studies of the Orion region, we explore possible explanations for the wide range in electron density 
derived from different ratios of trace and dominant species in the strongest neutral components, 
and we compare the elemental abundances and depletions with those found for other lines of sight. 
In an appendix, we briefly summarize some recent results on atomic data relevant to this work. 

2. Observations and Data Reduction 

2.1. 23 Ori Line of Sight 

The Bl V star 23 Ori [V = 5.00; E(B - V) = 0.11; vsaxi = 295 km s" 1 ] is located at (I, b) = 
(199?2, — 17?9), about 3 degrees south of 7 Ori (the northwest shoulder of Orion). The parallax 
3.39±0.87 mas measured by Hipparcos implies a distance of 295^qq 2 pc (ESA 1997), consistent 
with, but perhaps somewhat smaller than the value 430±100 pc inferred from the spectral type and 
apparent magnitude (using the absolute magnitudes of Blaauw 1963). Using Hipparcos parallaxes 
for a number of Orion stars, Brown, Walter, & Blaauw (1998) have suggested that Orion OB 
association subgroup la (of which 23 Ori is a member) lies at an average distance of about 330 
pc, roughly 100 pc closer than subgroups lb and lc. Morrell & Levato (1991) determined a stellar 
radial velocity v rac j = 37±6 km s _1 , somewhat larger than the previous estimate of 18 km s _1 due 
to Campbell & Moore (1928), and concluded that t> ra d does not vary. The line of sight to 23 Ori 
has a total neutral hydrogen column density log[./V(H I + 2H2)] = 20.74 cm -2 (Bohlin, Savage, & 
Drake 1978), with a relatively modest molecular hydrogen column density log[iV(H2)] ~ 18.3 cm" 2 
(Savage et al. 1977; Frisch & Jura 1980; this paper). 23 Ori may lie behind a detached wisp of Ha 
emission probably associated with Barnard's Loop; it does lie behind the high-velocity ionized gas 
dubbed "Orion's Cloak" by Cowie et al. (1979). While 23 Ori has been observed in connection 
with several surveys of optical and UV interstellar absorption- lines (e.g., Hobbs 1978b; Chaffee & 
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White 1982; Jenkins, Savage, & Spitzer 1986; Crane, Lambert, & Sheffer 1995), the line of sight 
has not previously been studied in detail. 

2.2. Optical Spectra 

High resolution, relatively high S/N observations of the optical absorption lines of Na I, Ca I, 
Ca II, K I, and several molecular species have been obtained using the 0.9 m coude feed telescope 
and coude spectrograph at Kitt Peak National Observatory, the 3.9 m Anglo-Australian Telescope 
and Ultra-High-Resolution Facility (UHRF) at the Anglo-Australian Observatory, and the 2.7 m 
telescope and coude spectrograph at the McDonald Observatory (Table ||). The references noted 
in Table || give details of the various instrumental setups and data-reduction procedures used to 
obtain and process the optical spectra. 

2.3. UV Spectra 

We observed 23 Ori with the HST GHRS during Cycle 1 (on 1992 April 14-15 and 1992 
September 23), as part of GO programs 2251 and 3993 aimed at determining the properties of 
individual interstellar clouds. Additional observations were obtained in Cycle 5 (on 1996 January 
13), as part of GO program 5896 studying high- velocity ionized gas. All observations reported here 
were obtained through the small science aperture, in order both to achieve the highest possible 
spectral resolution and to avoid complications due to the extended wings of the large science 
aperture instrumental function (for the pre-COSTAR spectra). Thirteen wavelength regions 
between 1740 and 2805 A, each covering 9-15 A, were observed using echelle-B at resolutions 
of 3.3-3.7 km s _1 (Table §). Owing to the (temporary) failure of side 1, and the consequent 
unavailability of echelle-A for obtaining high-resolution spectra at shorter wavelengths, ten 
wavelength regions between 1130 and 1680 A, each covering approximately 37-38 A, were observed 
with grating G160M at intermediate resolutions of 13-20 km s _1 (Table ^). The two ECH-A 
exposures, near 1260 and 1334 A, were obtained in Cycle 5. In order to obtain proper sampling 
and to reduce the effects of photocathode sensitivity variations, all exposures were divided into a 
number of slightly shifted subexposures. For the shorter wavelength G160M spectra, FP-SPLIT 
= FOUR was used to obtain four subexposures, successively shifted by typically 0.7 A; FP-SPLIT 
= DSFOUR was used for the ECH-A observations. For the longer wavelength ECH-B spectra, 
the WSCAN procedure was used to obtain three or four subexposures, each offset by 0.7-1.0 A. 
(See, e.g., Duncan & Ebbets 1989 for discussions of these observing strategies.) Additional G160M 
spectra of 23 Ori, obtained from the HST archive (GO program 2584), were used either to measure 
absorption lines not covered in our spectra or to increase the S/Ns for weak lines. 

Among the subexposures for each spectral region obtained in Cycle 1, the generally small 
relative offsets in the nominal wavelength scales supplied by STScI were determined by fitting the 
profiles of interstellar and/or stellar lines in the spectra. These relative offsets were typically less 
than 0.5 km s _1 for the ECH-B data, and usually less than about 10 km s _1 for the G160M data. 
The individual subexposures were shifted according to those relative offsets, sine interpolated 
(cf. Spitzer & Fitzpatrick 1993) onto a roughly twice-oversampled standard wavelength grid, and 
summed. Obvious deviant pixels in the individual subexposures (due to detector pattern features, 
bad pixels, radiation hits, etc.) were given zero weight in the sums; no other explicit detector 
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pattern noise identification and removal was performed, however. Adjacent pixels were then 
summed to yield approximately optimal sampling (i.e., two points per resolution element) in the 
final spectra. For the Cycle 5 ECH-A spectra, the STSDAS routines POFFSETS, DOPOFF, and 
SPECALIGN were used to determine the detector pattern and to align and sum the subexposures. 

Portions of the spectra near interstellar absorption lines of interest (usually including 150-250 
points) were normalized by fitting the continuum regions with Legendre polynomials of orders 
typically 7-15. For the Cycle 1 spectra, standard values for the on-order scattered light background 
(Cardelli, Ebbets, & Savage 1993) were subtracted from the normalized ECH-B spectra. The 
applicability of those standard values was confirmed via the strongly saturated Mg II lines at 2796 
and 2803 A, where the resulting relative intensity in the saturated line cores was within 0.1% of 
zero. The default (constant) backgrounds subtracted from the ECH-A spectra obtained in Cycle 5, 
however, left some residual flux (0.5-2.5%) in the otherwise clearly saturated cores of several lines; 
that residual flux was subtracted from those spectra. Except for the shortest wavelength G160M 
exposure centered near 1133 A, the continuum signal-to-noise ratios achieved for the G160M 
data are typically between 90 and 120, yielding 2a equivalent width upper limits or errors of 1-2 
mA for weak, unresolved absorption lines. The corresponding continuum S/Ns for the echelle 
data range from 30-40 for the initial set (I) of wavelength regions, from 60-100 for the second and 
third sets (II and III) of regions, and is about 120 for the region near the C II] line at 2325 A. The 
2a equivalent width errors or upper limits for the echelle data range from 0.6 mA near the C II] 
line to 1-2 mA for other weak lines. The quoted equivalent width errors include contributions 
from both photon noise and continuum placement (Jenkins et al. 1973; Sembach & Savage 1992), 
added in quadrature. 

The equivalent widths of absorption lines due to various atomic and ionic species found in 
the optical and UV spectra are listed in Table ||. For most of the lines, the equivalent widths 
correspond to the strong low- velocity absorption seen at km s _1 <J v <^ +30 km s ; absorption 
from intermediate- velocity gas at —45 km s _1 <J v ^ km s , blended with the low- velocity 
absorption, is also included where it is detected^} Equivalent widths and limits for the 
high- velocity (HV) components found at —108 km s _1 <J v & <^ —83 km s" 1 are noted separately in 
the last column (updating values given by Trapero et al. 1996). Lines observed at high resolution 
(FWHM <^ 4 km s _1 ) are noted by an "H" following the wavelength; dominant contributors 
to blended lines are marked with a "B" following the equivalent width. Minor contributors to 
blended lines are included in the table, but are referred to the dominant contributors for the total 
equivalent widths. The letters "P" and "S" denote primary and secondary transitions used in 
determining the column densities for each species. Equivalent widths and column densities for 
molecular absorption lines (including Copernicus data for H2) are given in Table ||. 

Offsets in velocity, relative to the zero point set by the high-resolution optical spectra of 
Ca II, are listed in the next-to-last column of Table ||. These offsets, determined in the analysis of 
the line profiles as described below, are consistent with those derived from the various wavelength 
calibration and "spybal" spectra obtained along with the stellar spectra. The offsets between the 
initial set of ECH-B spectra and the optical spectra are small, ranging from —0.1 to —1.7 km s , 
but the offsets for the second set of ECH-B spectra are significantly larger, ranging from —5.4 to 



4 In this study, we will use heliocentric velocities; the zero point for velocities with respect to the Local Standard 
of Rest is at V0 ~ +16.4 km s - . 
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—8.7 km s~ 1 . The zero point offset for the Al III spectra could not be determined by comparisons 
with the optical lines, since Al III is likely to be distributed significantly differently from those 
species. Since the offsets for each observing session generally showed systematic trends with time, 
however, we estimated an offset of —4.1 km s _1 for the Al III spectrum from the average of the 
offsets found for the immediately preceding spybal spectrum (—2.7 km s _1 ) and the immediately 
following Fe II A2249 spectrum (—5.4 km s _1 ). The relative velocity offsets determined from fits 
to the HV components are very similar to those derived independently from fits to the stronger 
low- velocity absorption. The average velocity offsets found between the G160M spectra and the 
optical and ECH-B data ranged from —40 km s _1 for the shorter wavelength settings to —25 km 
s _1 for the longer wavelength settings. Systematic variations of several km s _1 , with more negative 
values at shorter wavelengths, were obtained for the offsets within each G160M spectrum. 



2.4. Profile Analysis 

The normalized absorption-line profiles of various neutral and singly ionized species seen in 
the high-resolution optical spectra are shown in Figure |l[ Selected UV echelle and lower resolution 
G160M profiles, covering a larger velocity range, are shown in Figures [2] and ||. While the lower 
resolution (FWHM ^1.5 km s -1 ) optical and GHRS echelle profiles seem to show a single strong 
component at v ~ +23 km s _1 , the higher resolution (FWHM ~ 0.3-1.3 km s" 1 ) Ca II, Ca I, 
and K I spectra reveal (at least) two components dominating that strong observed absorption. 
A number of weaker (i.e., lower total hydrogen column density) components can be seen in the 
profiles of the stronger lines of various singly ionized species at km s~ x v <J +30 km s _1 . 
Additional weak intermediate-velocity (IV) components at —20 km s _1 <^ v <^ —4 km s _1 are 
apparent in the strong Mg II A2796 and A2803 lines; still weaker IV components at —43 km s _1 <^ 
v <^ —32 km s _1 can be seen in the very strong C II A1334 and Si III A1206 lines. The high-velocity 
(HV) gas comprising Orion's Cloak appears at —108 km s _1 <^ v <^ —83 km s _1 in the strongest 
lines of C II, N II, Mg II, Al II, Al III, Si II, and Si III; four components can be discerned in the 
echelle profiles of Mg II, Si II, and C II. Weak, broad absorption due to more highly ionized species 
(S III, C IV, Si IV, N V) is centered between —5 km s _1 and +6 km s . A schematic diagram of 
the observed absorption is given in Figure ||. 

Since the primary goal of this investigation is the determination of abundances and physical 
conditions for individual interstellar clouds, we have used the method of profile fitting in an 
attempt to discern and characterize the individual components giving rise to the observed 
absorption profiles (e.g., Welty, Hobbs, & York 1991). In this analysis, we attempt to obtain 
column densities (N), Doppler line width parameters [b = (2kT/m + 2v\ ) 1 ^ 2 ] and heliocentric 
velocities (v) for as many components as are needed (at minimum) to reproduce the observed 
line profiles. Comparisons among the various high-resolution optical line profiles reveal both the 
underlying component structure (number of components and their relative velocities, constraints 
on fe-values and relative abundances) and also some general characteristics of the individual 
components which together can be used to understand and model the lower resolution UV line 
profiles. In the analysis of the line profiles detailed below (and in the subsequent interpretation 
of the derived component column densities), we will make frequent use of the equation of 
photoionization equilibrium, T(X ) n(X°) = a(X°,T) n e n(X + ), where T is the photoionization 
rate and a is the recombination rate. For example, where abundances for several trace neutral 
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species are available for a given component, ratios of that equation for two different elements yield 
estimates for the relative total elemental abundances for the two elements (York 1980; Snow 1984; 
Welty et al. 1991). We note, however, that other processes besides simple photoionization and 
radiative recombination may affect the ionization balance for some elements. 

Where we had multiple lines from the same species, it was generally possible to find a 
consistent set of component parameters for fitting all the lines, within the uncertainties in the line 
profiles and/or equivalent widths. This consistency suggests both that the derived component 
model adequately describes this line of sight and that the adopted atomic data are accurate. For 
atomic transitions, we use the wavelengths, oscillator strengths, and damping constants tabulated 
by Morton (1991), except where more recent work has provided more accurate values (as noted 
in Table || and discussed in the Appendix). References for the corresponding molecular data are 
given in Table |j. 

The detailed component structures for the predominantly neutral gas at low velocities and 
for the predominantly ionized gas at intermediate and high velocities, derived primarily from 



the high-resolution optical and UV profiles, are described in Sections ^lj and ^lj below. The 
individual components are marked above the spectra in Figures [l], and ||] and are numbered 
sequentially, with the letters H, I, W, and S appended to associate them with the high-velocity, 
intermediate-velocity, weak low-velocity, and strong low-velocity component groups defined above 
and below. The absorption from the higher ions Al III, S III, C IV, Si IV, N V seen at relatively 



low velocities is discussed in Sections 4.1 and i.4 



2.5. iV(H I) and Reference Abundances 

We have adopted the total line of sight H I column density log[iV(H I)] = 20.74 cm -2 
determined by Bohlin et al. (1978) from Copernicus spectra, instead of the lower value 20.54 
cm -2 obtained by Diplas & Savage (1994) from IUE spectra. Profile fits to archival GHRS G160M 
spectra encompassing the H I Lyman a line, using the component structure derived for Zn II, yield 
log[iV(H I)] ~ 20.7 cm~ 2 . In addition, the lower value obtained by Diplas & Savage would imply 
essentially solar (or slightly higher) relative abundances for the typically very mildly depleted 
elements N, S, and Zn, as well as higher than usual abundances for C, O, and P. 

For calculating depletions, we use the meteoritic (CI carbonaceous chondritic) reference 
abundances of Anders & Grevesse (1989) and Grevesse & Noels (1993) for all elements except C, 
N, and O (which are solar photospheric values), as summarized in Table ||. For most elements, 
the meteoritic and solar photospheric abundances agree to within ±0.06 dex; the exceptions are 
P (A m -A p =0.12), CI (-0.23), Mn (0.14), Ga (0.25), Ge (0.22), and Sn (0.14) (though the listed 
uncertainties for photospheric CI and Sn are both 0.3 dex). For comparison with the depletions 
found toward 23 Ori, the Galactic "warm" and "cold cloud" interstellar depletion patterns listed 
in Table [5] have been adapted and updated from Jenkins (1987); the "halo cloud" pattern has been 
adapted from Savage & Sembach (1996) and Fitzpatrick (1996b). Slight differences from similar 
compilations in our previous papers (Lauroesch et al. 1996; Welty et al. 1997, 1999a) are due to 
the meteoritic reference abundances adopted here, to some newly-determined /-values (see the 
Appendix), and to recent results for the abundances of several elements (C: Sofia et al. 1997; N: 
Meyer et al. 1997; O: Meyer et al. 1998; Kr: Cardelli & Meyer 1997). The warm cloud depletion 
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for Ca is estimated in this paper (§ |3.2.3 ). Depletions for the heavy elements Cu, Ga, Ge, and Sn 



are based on data in Hobbs et al. (1993), Cardelli (1994), and Lambert et al. (1998). 

3. Predominantly Neutral Gas at Low Velocities 

3.1. Neutral Gas Component Structure 

For the absorption at relatively low velocities (v ^ km s , or |t>LSR| ^ 20 km s -1 ), O I, 
N I, and various singly ionized species show strong absorption. The main absorption component(s) 
near v ~ 23 km s _1 show the strongest lines both of various trace neutral species and of the 
corresponding dominant first ions. The absorption lines from the trace neutral species are 
remarkably strong, given the relatively modest total iV(H I); for example, some are stronger than 
the corresponding lines observed toward £ Oph, which has log[A r (H tot )] ~ 21.15 cm~ 2 . C I, Na I, 
and perhaps Mg I are the only trace neutral species detected in the somewhat weaker absorption 
at km s _1 <^ v <^ 20 km s . As expected, the component structure found for trace neutral 
species seen in the UV is very similar to that of Na I, K I, and Ca I (Welty et al. 1991). The 
absorption-line profiles for many of the dominant first ions of typically depleted elements (e.g., 
Fe II, Si II, Cr II, Ni II) appear somewhat similar to the Ca II profile, both toward 23 Ori and for 
a number of other lines of sight for which we have GHRS data. The latter similarity is somewhat 
remarkable, as Ca II is probably a trace ionization stage for many of those components, but is 
likely attributable to the offsetting effects of the Ca ionization balance and the Ca depletion for the 
different types of neutral (H I) clouds (Welty et al. 1996). The component structure found for K I, 
Ca I, and Na I was therefore used to model the UV profiles of other trace neutral species, while 
the Ca II component structure was used to model the UV profiles of other (generally dominant) 
first ions of comparable strength]^. The low-velocity component structures derived for species 
observed at high resolution are given in Table ^. 

3.1.1. Optical Spectra of Na I, K I, Ca I, and Ca II 

We begin the profile analysis with the high-resolution K I and Ca II profiles, where the 
detailed component structure is most clearly seen (Figure [[]). Independent fits to the highest 
resolution Ca II and K I profiles indicate that the strong absorption near +23 km s _1 arises 
primarily from two components, separated by about 2.10-2.15 km s _1 and of comparable strength 
(within a factor of ~ 2); several weaker outlying components are present as well. We will refer to 
these components (16S-19S) as the "strong low-velocity" (SLV) components. Six relatively weak 
blueward components (10W-15W) at km s _1 <^ v <^ +18 km s , and two others (20W-21W) at 
v ~ +28 and +30 km s , which we will call "weak low velocity" (WLV) components, are required 
to fit the profiles of Ca II, Na I, and some of the UV lines. The fits to the coude feed and higher 
resolution UHRF Ca II spectra are quite consistent for the SLV components at +20 km s _1 <J v <J 
+26 km s ; the UHRF spectrum is somewhat noisy, and not as useful for confirming details in 



5 In fitting the UV lines, the GHRS instrumental profile was assumed to be adequately described by a Gaussian, 
with the FWHM appropriate for each wavelength region taken from Duncan & Ebbets (1990); see column 4 and 
footnote b of Table 0. 
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the WLV absorption, however. 

Because the SLV Na I components are strongly saturated, we fit the Na I and K I profiles 
together, for various "reasonable" values of the ratio iV(Na l)/N(K I) for the two strongest 
components (17S and 18S). That ratio is typically within a factor of 2-3 of 70 in the nearby ISM 
(Welty et al. 1991). The Na I and K I column densities and 6-values for components 16S and 19S, 
and the overall velocity shift between the Na I and K I profiles, were determined via iterative fits 
for each set of those initial conditions. The resulting fits were evaluated on the basis of goodness 
of fit to the observed Na I profile and of the reasonableness of iV(Na l)/iV(K I) for those two 
outlying SLV components. iV(Na T)/N(K I) ~ 70 for the two strongest SLV components gave 
quite acceptable fits and very similar column density ratios for the two outlying components; 
iV(Na l)/JV(K I) ratios smaller than about 35 or larger than about 100 for the two strongest 
components yielded noticeably poorer fits. For the reasonable assumption that Na I and K I 
coexist, 6(Na I)/6(K I) will be in the range 1.0-1.3 (fully turbulent to fully thermal broadening). 
Since the best fits to the profiles were obtained for 6(Na I) ~ 6(K I), and since the 6-values are 
significantly larger than would be obtained for thermal broadening at T ~ 100 K (as estimated 
in Section |3.3.1 ), the broadening is apparently dominated by turbulence. The adopted Na I SLV 



component structure could be checked via future high-resolution observations of the Na I A3303 
lines. The Na I WLV components were fitted using the relative velocities determined from fits to 
the Ca II profile. 

The slight disagreement in relative velocities for SLV components 16S and 19S seen in 
the initial independent Ca II and Na I fits, the apparently larger 6-values for Ca II for those 
components, and initial fits to some of the UV lines all suggested the presence of additional 
structure. We therefore re-fit the Ca II and Na I profiles, with the SLV components fixed at 
the same relative velocities found for Na I and K I and with two additional WLV components 
suggested by the UV lines. The resulting Ca II 6-values are still slightly larger than those found for 
the corresponding components in Na I and K I, suggesting that the Ca II may occupy a somewhat 
larger volume (as found for a larger sample of clouds seen in both Na I and Ca II by Welty et al. 
1996)0. The weak Ca I absorption was fitted with two components having the same separation 
and 6-values as the two strongest K I components. 



3.1.2. Modelling the UV Spectra 

For the two strongest SLV components (17S and 18S), which are not resolved in the UV 
spectra, the column density ratio found for Ca I (1.0el0/0.6el0 ~ 1.7) was assumed to typify 
the relative abundances of the trace neutral lines of heavily depleted elements (Ca I, Fe I), while 
the column density ratio found for K I (1.37ell/1.61ell ~ 0.85) was assumed to hold for the 
trace neutral lines of lightly depleted elements (C I, Na I, S I, K I). Taking the ratio of the 
iV(Ca l)/iV(Ca II) ratios, for the same two components, yields the ratio of electron densities 
(0.026/0.010 ~ 2.6) in the two components (assuming similar T/ct). The iV(Ca I)/iV(Ca II) ratios 
also suggest that Ca II is the dominant ionization state of Ca in both components. The column 
density ratio for Call in the two strong components (3.9ell/6.3ell ~ 0.6) was therefore assumed 



6 The component structure given in Table |fj| for Ca II differs slightly from that given in our Ca II survey (Welty et 
al. 1996), which was based on fits to the Ca II spectra alone. 
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to apply for the dominant ions of other heavily depleted elements. Dividing the K I ratio by the 
electron density ratio (0.85/2.6 ~ 0.33) then yields an estimate for the relative abundances of the 
dominant ions of lightly depleted elements. For those two strongest components, the abundance 
ratios (r = N^s/Niss) of trace and dominant ionization states of elements typically depleted less 
than Ca but more than K were assigned intermediate values (i.e., 0.85 ^ r <J 1.7 for trace species 
and 0.33 <^ r 0.6 for dominant species), while maintaining the electron density ratio of 2.6. 
These predicted abundance ratios depend on the assumption that the 7V(X I)/iV(X II) ratios for 
different elements X all depend in the same way on n e in the two strongest SLV components — as 
would be the case if photoionization and radiative recombination dominate the ionization balance. 
We will find in Sections |3.3.3 and 5.2, however, that different neutral/first ion ratios yield different 
values for n e (for a given T and radiation field) — suggesting that additional processes may be 
affecting the ionization balance. In particular, the n e inferred from the iV(Ca I)/iV(Ca II) ratio 
is often significantly higher than the values derived from other such ratios. The relative column 
densities for components 17S and 18S may therefore differ from the values assumed for the profile 
fits. We note, however, that consistent abundance ratios were obtained from fits in which Nns 
and iVi85 were allowed to vary - both for the lightly depleted Zn II and for the more severely 
depleted Fe II. Since at least the two strongest SLV components appear to consist of relatively 
cold, but turbulent gas, the b- values for most species are very similar for each component. For the 
apparently warmer WLV components, the 6-values for various species were derived from the WLV 
Ca II 6-values assuming a temperature of 3000 K (S 



3.1.3. UV Echelle Spectra 

For lines due to most trace neutral species and for the weaker lines due to some first ions, 
the WLV components were not detected in the UV echelle profiles, so that it is very difficult 
to unambiguously determine both the individual component column densities for the SLV 
components and the overall velocity offset, given the typical ~ 1 km s _1 inherent uncertainties 
in the GHRS echelle velocities. In such cases, the optical SLV column densities were uniformly 
scaled to fit the UV profiles, after constraining the ratio of the two strongest components as 
just described — in effect yielding only an overall SLV column density ratio and velocity offset. 
While this uniform scaling gave an acceptable fit in all cases, the UV line profiles do permit 
some component-to-component variations. Where the WLV components are present in the echelle 
spectra, however, they provide additional constraints on the overall velocity offset, and it was thus 
possible to estimate individual component column densities for Zn II (A2026, A2062), Fe II (A2344, 
A2374), Si II (A1808), and S II (A1259), as listed in Table ||. Separate simultaneous fits were made 
to the three Fe II lines and to the two Zn II lines observed with the echelle. In the S II A1259 
profile, the SLV components are saturated, and were thus poorly constrained in the profile fits. 
For "reasonable" values of SLV N(S II) and the velocity zero-point offset (constrained by Si II 
A1260 and C I A1260), however, N(S II) is well determined for all of the WLV components except 
15W and 20W, which are adjacent to the SLV components. Since the WLV N(S n)/JV(Zn II) 
ratios suggest that sulfur is undepleted in the WLV components, we estimated the total N(S II) 
by using the curve of growth defined by the other dominant, little-depleted species and by the 
derived Zn II component structure (§ |3,2.4 ). Integrated over the line of sight, sulfur appears to be 



essentially undepleted, within ±0.2 dex. The very strong C II A1334 and Mg II A2796 and A2803 
lines are saturated over the entire WLV-SLV velocity range, and exhibit relatively weak, but 
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definite damping wings (Figures || and ||). We obtained total WLV + SLV column densities for 
C II and Mg II by using the component structures determined for Zn II and Si II, respectively, to 
fit the damped profiles — uniformly scaling the WLV and SLV column densities until a "smooth" 
continuum was achieved (as often done in fitting damped lines of H I). For those fits, the velocity 
zero points were constrained by requiring alignment of the high-velocity components in the echelle 
spectra of C II, Mg II, and Si II, which show very similar HV component structure. The resulting 
WLV + SLV column densities for C II and Mg II are consistent with the limit determined from 
the weak C II] A2325 line and with the value obtained from the weak Mg II A1239 doublet, 
respectively, using the revised /-values for these weak transitions discussed in the Appendix. In 
view of the complexity of the component structure and the relatively low S/N obtained for some of 
the profiles, the derived individual component parameters listed in Table ^ must still be regarded 
as somewhat uncertain. 

3.1.4. G160M Spectra 

For the typically higher S/N but lower resolution G160M spectra, individual component 
parameters could not be reliably derived. We therefore used the component velocity structure and 
relative column densities derived from the optical and UV echelle spectra to model the G160M 
spectra. For the WLV and SLV gas, the echelle profiles of Zn II, Si II, and Fe II (representing 
a range of depletion behavior for species dominant in H I regions) were especially useful for 
interpreting the lines of other dominant neutral and singly ionized species observed with G160M. 
Where the WLV components could be seen (e.g., as a blueward wing in the absorption- line profiles 
in some of the longer wavelength G160M spectra), the total SLV and WLV column densities could 
be scaled separately, guided by fits to the echelle spectra of similarly depleted species, to fit the 
observed G160M profiles. Where the WLV components could not be separately discerned (e.g., 
for weak lines and for the shorter wavelength, lower resolution G160M spectra), relative SLV and 
WLV column density scalings from species with similar overall depletion were used, essentially to 
fit the observed equivalent width. In such cases, the uncertainties in N correspond to the range in 
values required to yield W\ =L la. The derived WLV and SLV column densities are typically less 
reliable for species (such as Al II) which have only a single strong line in the G160M spectra. 

3.2. Abundances and Depletions for SLV and WLV Gas 

In Table g, we list column densities or upper limits for the predominantly neutral (H I) gas 
- for the total line of sight and for the SLV and WLV component groups at km s _1 <J v <^ +30 
km s . The code in the last column of the table indicates how the component group column 
densities were determined, as discussed in the previous section. Species whose WLV and SLV 
column densities were obtained via scaling from similarly depleted species (typically observed at 
higher resolution) are given in parentheses in Table |?]. For several species (e.g., Mn I, Kr I), we 
have averaged spectra from several weak transitions in order to obtain either a detection or a more 
stringent limit. Mn I is detected (tentatively) for the first time in the diffuse Galactic ISM toward 
23 Ori. 

In Figure ||, we show several column density ratios for the individual WLV and SLV 
components to illustrate some of the differences between those two sets of neutral components. At 
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the right in some of the panels, we show the ratios for the meteoritic reference abundances and for 
the three representative Galactic interstellar cloud types listed in Table B The iV(Na l)/iV(Ca II) 
ratio provides the clearest distinction between the WLV and SLV components, though the 
N(Si n)/JV(Zn II), N(Fe H)/iV(Zn II), and iV(Fe II)/JV(S II) ratios reveal some differences in 
depletion (§§ |3.2.1| and |5.3.3|) . As noted above, the ratio iV(Ca T)/N(C& II) suggests that Ca II is 
the dominant ionization state of Ca in the two strongest SLV components (17S and 18S) (Welty 
et al. 1996, 1999b). The smaller iV(Ca Il)/JV(Fe II) ratios for the outlying SLV components (16S 
and 19S), however, may indicate that Ca II is a trace species there. Similarly, comparing the ratio 
iV(Na l)/iV(Ca II) with iV(Na I) suggests that Ca II is likely a trace species in the weaker SLV 
and WLV components (Welty et al. 1996). The iV(Na l)/iV(Zn II) ratio provides constraints on 
the electron densities in the individual WLV components (if a typical interstellar radiation field is 
assumed). Finally, we note that the slightly sub-solar N(Zn ll)/N(S II) ratio found for the WLV 
components in the profile fits (and the essentially identical value obtained for the SLV components 
from the curve of growth) are consistent with the results of Fitzpatrick & Spitzer (1997) for the 
predominantly halo clouds toward HD 215733. 

While the Al III, Si III, and N II profiles indicate that there is apparently some ionized 
(i.e., H II) gas present at km s _1 <^ v <J +30 km s _1 (Figure ||), comparison of the column 
densities with those of Al II, Si II, and N I at similar velocities implies that the SLV and WLV 



gas is predominantly neutral. For example, the limit on iV(N II) (§§ 4.1 and 4.4) suggests that 
the low velocity gas is less than about 5% ionized (though the limit on the fractional ionization 
is somewhat higher for the WLV gas and somewhat lower for the SLV gas). The relative column 
densities of O I, N I, and C II, very similar to those typically found in H I gas, also suggest that 
the SLV and WLV gas is mostly neutral, since only C II would be present in any ionized gas. As 
we will show below (§ ||), the various IV and HV components at v <^ km s" 1 are apparently 
predominantly ionized. It thus seems reasonable to ascribe the measured N(H I + 2H2) to the 
SLV and WLV gas. If the contribution of any H II gas to the observed Zn II and S II absorption 
is negligible, then Zn and S are nearly undepleted in the SLV and WLV gas (Dz n = —0.12 dex; 
Ds ~ 0.0 dex). We may then use the Zn II and S II column densities found for the individual SLV 
and WLV components to infer neutral hydrogen column densities for those individual components 
(assuming the same small level of depletion in each). Unsurprisingly, the SLV components contain 
the bulk of the neutral gas, with total log[A r (H)] = 20.71 cm -2 . The WLV components have a 
total log[JV(H)] = 19.61 cm~ 2 . 

The depletions of various elements listed in Table ^ and plotted in Figure |6] for the SLV and 
WLV components (integrated over four and seven components, respectively) are with respect to 
the neutral hydrogen column densities inferred from Zn II and S II. In Figure [f| the elements are 
ordered according to the condensation temperature, as tabulated by Wasson (1985). The average 
depletions for the line of sight are very close to those found for the SLV components, which contain 
most of the total H I. Where we do not have independent values for the WLV and SLV depletions, 
we have therefore plotted only the average line-of-sight values (with an "s"). For comparison in 
Figure ||, the Galactic "warm" and "cold cloud" depletion patterns (Table ||) are plotted as dotted 
and solid lines, respectively. The £ Oph "component A" (warm, low density gas near vq ~ —28 
km s _1 ) and "component B" (cold, dense gas near vq ~ —14 km s^ 1 ) depletions are derived from 
data in Savage, Cardelli, & Sofia (1992), Hobbs et al. (1993), Cardelli et al. (1994), Sofia et al. 
(1994), Lambert et al. (1998), and from our own analyses of various optical and archival GHRS 
spectra. We have adopted log[A r (H tot )] = 19.95 cm -2 for £ Oph component A by assuming D(Zn) 
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~ —0.2 dex — as found for warm, low density gas by Sembach et al. (1995). 

3.2.1. SLV Depletions 

For most elements, the 23 Ori SLV depletions (Table [7| and Figure |6|) appear to be less severe, 
by factors of 2-4, than the typical cold, dense cloud depletions or the depletions found for the 
strong component B blend seen toward Q Oph. It is conceivable that such an intermediate depletion 
could arise from a mixture of several components, each characterized by either "typical" warm 
or cold cloud depletions. As noted above, the depletions in the two strongest SLV components 
(17S, 18S) differ by factors 2 for different elements, and the weaker SLV components (16S, 
19S) may be warmer and less dense than the two strong SLV components. Where we have 
echelle data, however, the line profiles do not seem to permit a mixture extreme enough to 
yield the overall factor of 2-4 (see Fig. ||). While the iV(Fe II)/iV(Zn II) ratios suggest that the 
depletions for component 16S may be somewhat less severe than those for components 17S-19S, 
the N(N& l)/./V(Ca II) ratio is somewhat higher for both the outlying SLV components (16S and 
19S) than for the WLV components. The possible presence of a small amount of ionized gas would 
affect the depletions by less than 10% (0.04 dex). Interestingly, the SLV depletions are much closer 
to the cold cloud pattern when viewed as [X/Zn] = (X/Zn)^ — (X/Zii)q, since Zn also appears 
to be somewhat less severely depleted than in the cold cloud pattern. 

3.2.2. SLV Molecules 

Absorption from the molecular species H2, CH, and CH + is also present at velocities 
corresponding to the SLV components. The H2 column densities for rotational levels J=0 through 
J=5 listed in Table ||] have been derived by comparing the observed {Copernicus) H2 line profiles 
with model profiles based on the adopted SLV component structure for Na I — using either all 
four SLV components or just the two strongest ones. While Frisch &: Jura (1980) used different 
6-values for the odd and even rotational levels to obtain estimates for the H2 column densities, 
our fits do not suggest any obvious systematic difference in component structure between the odd 
and even levels (though the resolution and S/Ns are rather poor). The column densities in the 
lowest four rotational levels yield an ortho- to para-H_2 ratio of ~3, similar to the values found 
by Spitzer, Cochran, & Hirshfeld (1974) for lines of sight with somewhat lower iV(H2). The CH 
absorption in the SLV components toward 23 Ori is remarkably strong — a factor of at least 40 
stronger than would be expected from the usually fairly tight correlation between iV(CH) and 
iV(H 2 ) (e.g., Danks, Federman, & Lambert 1984)[|. The 7V(CH)/iV(CO) ratio is also higher 
than in most other lines of sight. Federman, Welty, & Cardelli (1997) have proposed that the 
CH observed toward 23 Ori is due to the same (poorly understood) non-equilibrium process(es) 
responsible for the observed CH + , rather than to the equilibrium gas-phase chemical reactions 
responsible for the diatomic molecules (e.g., CH, CN, C2) observed in many other diffuse cloud 
lines of sight. While the absorption from both CH and CH + may be acceptably fitted with a single 
component (Table ||), the presence of two strong SLV components in the atomic species and the 
possible asymmetry in the CH + profile (Figure |l]) suggest that two components may be present in 



7 The even higher JV(CH) reported by Crane et al. (1995) appears to be in error (P. Crane, private communication). 
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CH and CH + as well (see also Crane et al. 1995). The relatively broad components required even 
for two-component fits and the lambda doubling present for CH make more complex fits not well 
constrained, however. 

3.2.3. WLV Depletions 

The 23 Ori WLV depletions (Table @ and Figure g) are somewhat less well known for the 
species whose relative SLV and WLV column densities were constrained to be the same as for Zn II, 
Si II, or Fe II, based on similarities in average line-of-sight depletion. Overall, however, the WLV 
depletions seem to be less severe by an additional factor of 2-3 than the SLV depletions, and are 
similar to the typical warm cloud and £ Oph component A values. We obtained the Ca depletion 
for the WLV gas by using the N(C l)/N(C II), iV(Na I)/JV(Na II), and iV(Mg l)/iV(Mg II) 
ratios to estimate n e , then using n e to derive the column density of the presumably dominant 
Ca III. For temperatures from 100 K to 3000 K, n e ranges over roughly an order of magnitude 
(§ |3.4p , but D(Ca) ~ —2.4 dex throughout — roughly a factor 8 more severe than that of Fe 
in the WLV gas. A similar analysis applied to the £ Oph A components yields very similar 
values for the Ca depletion there, and Ca also seems more severely depleted than Fe in the 
component groups exhibiting "warm cloud" depletions toward SN 1987A (Welty et al. 1999a). If 
the iV(Ca Il)/iV(Ca III) ratio depends on n e in the same way as do the other three ratios, then 
Ca may generally be more severely depleted than Fe in warm, diffuse gas, contrary to what is 
often assumed. Figure 4 in Jenkins (1987), for example, shows D(Ca) ~ D(Fe) ~ —1.2 dex for 
log(< nu >) = —1-5 cm~ 3 , where the Ca depletion is based largely on relatively low resolution 
optical and IUE spectra obtained by Phillips, Pettini, & Gondhalekar (1984). Improved atomic 
data, allowance for the (very mild) depletion of Zn, and better data for several stars (especially a 
Vir; Hobbs 1978a) make the Ca depletion at that low < «h > more severe by about a factor of 5. 
We have therefore adopted a "warm cloud" Ca depletion of —2.0 dex in Table || and in Figure |6|. 

3.2.4- Empirical Curves of Growth 

Given the column densities for various neutral and singly ionized species in the SLV and WLV 
gas derived from fits to the observed line profiles and/or equivalent widths, we may construct an 
empirical curve of growth for the neutral gas along the line of sight (Figure 0) . Curves of growth 
have often been used to estimate line-of-sight column densities from lower resolution spectra, 
usually by comparing the observed log(W>/A) vs. log(/A) with theoretical curves derived for 
assumed "representative" 6-values. In the figure, we have used separate symbols for species that 
are trace ions in the H I gas (squares), species that are dominant ions of mildly depleted elements 
(e.g., C II, O I, Zn II) (triangles), and species that are dominant ions of more severely depleted 
elements (e.g., Si II, Ti II, Fe II) (circles). The solid lines are the curves that would be expected 
from the adopted component structures for Na I (trace species), Zn II (dominant, mildly depleted), 
and Fe II (dominant, severely depleted); the diagonal line shows the relationship for unsaturated 
lines. It is apparent that the three groups define distinct curves of growth, reflecting the differences 
in distribution and physical properties among the various SLV and WLV components. The curve 
defined by the trace neutral species, which are concentrated in the two strongest SLV components, 
is the first to show signs of saturation in the absorption lines as log(Nf\) increases, and lies 
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below the other two curves for 7.5 cm -1 ^ log(iV/A) <^ 10.5 cm -1 . The curve occupied by the 
dominant, depleted species, which are less severely depleted (and thus relatively stronger) in 
the WLV components, lies above the other two curves — reflecting the effectively somewhat 
broader distribution of the dominant, depleted species. The curve defined by the dominant, mildly 
depleted species falls in between the other two curves. Note that the strong lines of N I and O I, 
whose column densities were derived from weak semi-forbidden transitions, fall nicely along the 
predicted curve for dominant, mildly depleted species [near log(NfX) ~ 9.8-10.6 cm' 1 ], even 
though the curve is based on the component structure found for the significantly weaker Zn II 
lines. The three triangles nearest log(A r /A) ~ 9 cm' 1 represent the S II A1250, A1253, and A1259 
lines. The strongest lines of Mg II and C II exhibit both damping wings and weak IV absorption, 
and so lie somewhat above the predicted curves based on Fe II and Zn II, respectively. 

While differences in the curves defined by trace and dominant species have been noted before 
(e.g., Morton 1975; Spitzer & Jenkins 1975), dominant species with significantly different levels of 
depletion have often been assumed to occupy the same curve. Toward 23 Ori, that assumption 
could lead to errors of factors of 2-3 in inferred column density for some species, depending on 
the particular absorption lines available, even though the SLV and WLV depletions differ by 
only about 0.4 dex for most typically depleted elements. Furthermore, the curves of growth for 
both trace neutral species and dominant, mildly depleted species exhibit "kinks", at log(iV/A) ~ 
9.5 and ~ 8.5 cm' 1 , respectively — i.e., neither curve can be characterized by a single effective 
6-value. The potential presence of "kinks" in curves of growth complicates the use of such curves 
for determining "mutually consistent" /-values for lines of a particular species. 



3.3. Physical Conditions in SLV Gas 

3.3.1. SLV Temperature 

The strong absorption due to trace neutral species, the presence of moderate amounts of 
H2, fairly large values for A^(Na l)/./V(Ca II) ~ 10-30, and indications that Ca II is the dominant 
ionization state of Ca all suggest that the two strongest SLV components (17S and 18S, near +23 
km s _1 ) are relatively cold and dense. The apparent equality in b for K I and Na I, which differ 
by a factor 1.7 in atomic weight, also suggests that these components are due to relatively cold 
gas, with the absorption line widths dominated by turbulent broadening or bulk motions. The K I 
6-values for the two strongest SLV components yield strict upper limits for T of about 1870 K and 
4200 K, respectively, though the temperatures are actually much smaller if turbulence dominates 
the broadening. The slightly larger 6(Ca II) may indicate that Ca II occupies a larger volume, 
characterized by larger temperature and/or turbulence (Welty et al. 1996). 

For IJ_2 column densities greater than about 10 17 cm' 2 , the relative populations of the J=0 
and J=l rotational levels of H2 are thermalized via collisions with protons, and we can estimate 
the local kinetic temperature T from the ratio N(1)/N(0): 

^ = ^ exp(-£ iAT 01 ) = 9 exp(-170 K/T 01 ). 
N{0) go 

The temperature inferred from iV(l)/iV(0) lies in the range 65-150 K, depending on the exact 
component structure assumed for the H2 (cf. Frisch & Jura 1980); we will adopt T = 100 K for 
the two strongest SLV components. The analogous "temperature" inferred from the N(A)/N(2) 
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ratio is about 350 K — similar to values found by Spitzer et al. (1974) for other lines of sight and 
consistent with the relative level populations predicted by Spitzer & Zweibel (1974) for ran = 10 
cm -3 , T = 80 K, and an "average" interstellar radiation field. 



3.3.2. SLV Thermal Pressure and Local Density 

The relative populations of the fine structure levels in the ground electronic term of C I 
can yield information concerning the product n^T (i.e., the thermal pressure) in the SLV clouds 
(Jenkins & Shaya 1979; Jenkins, Jura, & Loewenstein 1983). If we have independent constraints 
on the temperature, then we can determine the local hydrogen density. Taking N(C I) = 8.5±1.0 
x 10 14 cm" 2 , N(C I*) = 9.4±1.8 x 10 13 cm" 2 , and N(C I**) = 1.4±0.2 x 10 13 cm" 2 , we have 
N(C I**)/JV(C l)tot = 0.015±0.003 and N(C l*)/N(C l) to t = 0.10±0.02. Those ratios imply 
log(n H T) ~ 3.1±0.1 cm~ 3 K for T ~ 100 K and the WJ1 radiation field — at the low end of the 
range of values reported by Jenkins et al. (1983). Inclusion of a "turbulent" pressure term, using 
the vt obtained from the K I 6-values for the two strongest SLV components, would increase the 
pressure by no more than a factor of about 2. If T ~ 100 K, as derived from H2, then tih ~ 10-15 
cm -3 . The total thickness implied for the two main clouds is thus about 12-16 pc — comparable 
to the apparent thickness of the 21 cm emission filaments, at SLV velocities near 23 Ori, which 
define the shell surrounding the Orion-Eridanus bubble in Brown et al. (1995; see discussion in 



Section 5T). A relatively long path length through the expanding shell might also explain why 
the 6-values for SLV components 17S and 18S are significantly larger than the thermal widths 
at T ~ 100 K — i.e., the apparent dominance of "turbulent " broadening could be due to slight 
differences in velocity through the shell. 



3.3.3. SLV Electron Density and Fractional Ionization 

If we assume a "typical" interstellar radiation field (e.g., the WJ1 field) and T ~ 100 K, we 
can use the equation of ionization equilibrium, T(X°) n(X°) = a(X°,T) n e n(X + ), to compute 
local electron densities n e for various pairs of neutral and singly ionized species (Table |8|). We 
have assumed depletions of —0.5 dex for Na and K (slightly less severe than the cold cloud values 
in Table ||) to estimate SLV column densities for Na II and K II. The n e estimated from Na and 
K are then consistent with those obtained from Mg and Si, whose ionization balance should be 
dominated by simple photoionization and radiative recombination in the SLV gas (§ |5.2j ). All 
other SLV column densities in Table || are measured values. The photoionization rate V depends 
on the ambient radiation field and the optical depth (e.g., Ay) within the cloud, and can be 
parameterized as F = C exp(— a Ay + (3 Ay), where C is the unattenuated rate and a and 
(5 (typically in the range 3 to 5) are fit coefficients derived for a chosen grain model (e.g., van 
Dishoeck 1988). Toward 23 Ori, the small total E(B — V) = 0.11 suggests that the attenuation 
will generally be small within the individual clouds along the line of sight — probably not by more 
than 25±5% (for different species) at the centers of the two strongest SLV clouds. We have listed 
the unattenuated rates for the WJ1 and Draine (1978) radiation fields in Table ||| The radiative 
recombination coefficient a depends on the temperature, and can be approximated by a = a Q 
(T/T )~ v for T ^ 10 4 K, where r\ ~ 0.6-0.7 for most species of interest (Pequignot & Aldrovandi 
1986). While the true values for n e will depend on the actual temperature and radiation field, the 



-18- 



values listed in Table || (for T = 100 K and the WJ1 field) can be compared in a relative sense, 
due to the similarities in the response of different species to changes in those parameters. [These 
n e are average values for the two strongest SLV components; recall that ra e (17S) ~ 2.6 n e (18S), as 
estimated from N(Ca l)/iV(Ca II).] We find that the computed average values of n e range over a 
factor of about 25 — from about 0.04-0.05 cm -3 for S, Mn, and Fe to about 0.09-0.26 cm -3 for 
C, Na, Mg, Si, and K to about 0.95 cm -3 for Ca. We will adopt for the present n e ~ 0.15±0.05 
cm -3 , from Na, Mg, Si, and K. With the ran derived above, the fractional ionization n e /nu ~ 0.01 
- significantly larger than the value expected if most electrons come from photoionization of C 
(i.e., ~ 2 x 10~ 4 ) — which suggests that hydrogen is partially ionized (at a level ~ 1%) in the 
SLV gas. We will discuss both the large range in n e and the apparently high fractional ionization 



in the SLV gas further below (§ |5.2[ ). As discussed in Section 5.1, it is unlikely that the ambient 
radiation field is weaker than the "typical" galactic field (in which case the true n e and n e /nn 
would be lower than the values computed here). 

We may use the n e inferred from photoionization equilibrium, together with the T and nn 
derived above, to predict the column densities of the excited fine-structure levels of C II and Si II 
in the predominantly neutral SLV gas. The equation of excitation equilibrium, for element X = C 
or Si, is (Bahcall & Wolf 1968; York & Kinahan 1979; Fitzpatrick & Spitzer 1997): 

N(X IP) _ n e7l2 (e)+n Hll2 (H) 



N(X II) A 21 + n el2 i{e) + nnl2i{Hy 

Here the collisional excitation rate for electrons is 

, . 8.63 x 10~ 6 n Blo/fcT , i 
712(e) = /2 fi 12 e ~ El2 / kT cm 3 s x . 

For C II, gi = 2, £12 = 0.0079 eV, the collision strength Q 12 ranges from 1.80 at T = 100 K 
to 2.26 at T = 3000 K to 2.85 at T = 8000 K (Hayes & Nussbaumer 1984), the collision rate 
for H atoms 712(H) ranges from about 5 x 10~ 10 cm 3 s -1 at T = 100 K to about 1.2 x 10~ 9 
cm 3 s _1 at T = 3000 K (York & Kinahan 1979), and the radiative decay rate A 21 = 2.29 x 
10~ 6 s- 1 (Mendoza 1983). For Si II, g x = 2, E 12 = 0.036 eV, n 12 ~ 5.66 (within ± 0.12 for T 
between 100 K and 20,000 K; Keenan et al. 1985), 712(H) ranges from about 2 x 10~ n cm 3 s _1 
at T = 100 K to about 10~ 9 cm 3 s" 1 at T = 3000 K (York & Kinahan 1979), and A 21 = 2.17 x 
10 s _1 (Nussbaumer 1977). For n e = 0.15 cm -3 , ran = 12 cm~ 3 , and T = 100 K, the excited 
state populations are dominated by electron collisional excitation and radiative decay. Using the 
SLV column densities of C II and Si II listed in Table |7|, we would predict log[iV(C II*)] ~ 15.3 
cm -2 and log[iV(Si II*)] ~ 10.8 cm -2 . The best fit to the observed C II* A1335 profile, using the 
component structure found for S II and Zn II, yields log[iV(C II*)] ~ 14.2 cm -2 , though values up 
to about 15.0 cm~ 2 may be possible if some residual background remains at SLV velocities. The 
predicted SLV column density of C II* is thus a factor of 2 to 13 larger than the observed column 
density. The predicted SLV log[iV(Si II*)] is slightly greater than the upper limit (10.7 cm" 2 ) 
obtained from scaling the derived Si II component structure to fit the weak A1264 absorption. If we 
reduce n e by a factor 10 (a possibility discussed in Section ^ below), the predicted SLV column 
densities become log[A^(C II*)] ~ 14.6 cm -2 and log[iV(Si II*)] ~ 9.9 cm -2 , more consistent with 
the observed values. Similarly, Spitzer & Fitzpatrick (1993) found the n e inferred from C II* 
fine structure excitation to be systematically lower than those determined from Na ionization 
equilibrium in several components toward HD 93521. The C II* observed at SLV velocities thus 
may arise in the predominantly cool, neutral (H at most ~ 1% ionized) SLV gas. 
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The physical conditions derived for the SLV gas (averaged over the four SLV components) 
are summarized in Table pi Note that the analysis of the C I fine structure levels (from which nn 
was obtained) assumes that the relative column densities for all three C I fine structure levels in 
components 16S-19S are the same as those found for K I and that each of T and nn are equal in 
those clouds. The fractional ionization would thus be different for components 17S and 18S by 
the factor 2.6 difference in n e (§ |3.1.2j ). To obtain equal n e /nn, one could, in principle, adjust the 
relative column densities for C I* and C I**, so that instead T and/or nn would be different for 
those two components. The required column density adjustments would yield velocity offsets for 
C I* and C I** that differ by about 0.4 km s _1 from that of C I for the A1260 multiplet observed 
with ECH-A. 



3.4. Physical Conditions in WLV Gas 

The weaker absorption from most trace neutral species (C I, Na I, and perhaps Mg I are the 
only such species detected), the smaller values for iV(Na T)/N(C& II) ~ 0.3-2.8, the generally 
somewhat milder depletions, and the larger 6-values for the individual Ca II components suggest 
that the WLV gas is warmer and/or of lower density than the SLV gas. We assume that the WLV 
gas is predominantly neutral, as fits to the G160M profiles of the strong O I A1302 and N I A1200 
lines seem to indicate that the WLV column densities of those two species are unlikely to be much 
smaller than the values listed in Table 0, though uncertainties as to the exact velocity offsets make 
it difficult to determine strong constraints. The ratio N(Mg l)/JV(Na I) ~ 5 (larger by less than 
a factor 2 than for the SLV gas) suggests that T <J 5000 K, since Mg I would be significantly 
enhanced at higher temperatures due to dielectronic recombination^. In the analysis below, we 
will consider possible WLV temperatures of 100 K and 3000 K, representative of cold and warm 
gas, respectively. 

Analysis of the WLV C I absorption yields relative fine structure level populations 
JV(C I**)/N(C l)tot = 0.16±0.04 and JV(C I*)/N(C l) to t = 0.40±0.08. While the nominal ratios 
do not fall along the curves in Jenkins & Shaya (1979) for uniform clouds — which would suggest 
that the individual WLV components are characterized by different physical properties — the 
region permitted by the uncertainties does substantially intersect the curves. If the WLV gas is 
cold (T ~ 100 K), then log(n H T) ~ 4.3 cm~ 3 K and n H ~ 200 cm" 3 . For warmer gas (T ~ 3000 
K), log(nuT) ~ 4.7-4.8 cm _3 K and nu ~ 15-20 cm -3 . Estimates of the average WLV n e from 
the ratios of N(C T)/N(C II), iV(Na l)/JV(Na II), and N(Mg l)/N(Mg II) (assuming a depletion 
of -0.3 dex for Na) range from about 0.02 cm~ 3 for T ~ 100 K to about 0.2 cm" 3 for T ~ 3000 
K (Table ||; bottom). For T = 100 K, we thus would have n e /ni{ ~ 0.0001 — i.e., lower than 
the minimum value expected from photoionization of C alone. For T = 3000 K, we would have 
n e /nn ~ 0.01 — comparable to the fractional ionization inferred for the SLV components. Since 
the higher temperature also yields better agreement with the observed C II* and Si II* absorption 
(below), we will adopt T ~ 3000 K, but note that the gas could be somewhat cooler. For T 
~ 3000 K, the density of the WLV gas would thus be slightly greater than that of the colder, 
higher column density SLV gas, the thermal pressure would be significantly higher, and the total 



8 The WLV N(Mg I) is based on the possible weak absorption in the Mg I A2026 profile in the same velocity range 
as the strongest WLV Na I components, and is not likely to be significantly underestimated. 
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thickness (0.7-0.9 pc) would be much smaller. If the temperatures and radiation fields are similar 
among the WLV components, the ratios iV(Na l)/iV(Zn II) suggest that n e may be higher by 
about a factor 3 for components 11W and 12W than for the other WLV gas (Fig. ||). 

As for the SLV gas, we may compare the observed C II* and Si II* absorption at WLV 
velocities with predictions based on the WLV n e inferred for different T. Fits to the ECH-A 
C II* A1335 profile, using the WLV and SLV component structure found for S II and Zn II, yield 
log[JV(C II*)] ~ 14.3-14.6 cm" 2 for the WLV gasg. Using the WLV N{C II) in Table % for T 
= 3000 K, n e = 0.2 cm" 3 , and n H = 20 cnr 3 , we predict log[7V(C II*)] ~ 14.25 cm" 2 , while 
for T = 100 K, n e = 0.02 cm" 3 , and n H = 200 cm" 3 we predict log[iV(C II*)] ~ 14.5 cm" 2 . 
For the same two combinations of T, n e , and nn, we predict log[iV(Si II*)] ~ 11.3 cm~ 2 and 
10.0 cm -2 , respectively, compared to the observed log[iV(Si II*)] ~ 11.5 (most of the observed 
Si II* absorption appears to be at WLV velocities). If the gas is cold, the dominant processes 
are excitation by H atoms (because of the high nji) and radiative decay; collisions with electrons 
dominate the excitation in warmer gas. These comparisons indicate that the C II* and Si II* 
absorption at WLV velocities could arise in the predominantly neutral WLV gas, if that gas is 
warm (T ~ 500-5000 K). While N(C H*)/JV(C II) is of order 0.01 for most of the SLV and WLV 
components [assuming N(C Il)/iV(S II) ~ 10], the ratio appears to be significantly higher for the 
WLV components 10W, 11W, and (perhaps) 12W — which may indicate a higher n e in those 
components, as suggested also by the A^(Na l)/iV(Zn II) ratio (above). 

The physical conditions derived for the WLV gas (averaged over the seven WLV components, 
for the two representative temperatures) are summarized in Table ||. Much of the characterization 
of the WLV components depends on the C I fine-structure level populations, which were derived 
from fits to the apparent blueward wing in the G160M profiles of the A1560 and A1656 multiplets. 
The WLV column densities are consistent with the limits obtained from the weaker A 1260 multiplet 
observed with ECH-A, but higher resolution data for one of the stronger multiplets are needed 
to confirm the relatively high C I excited state populations in the WLV gas. We note, however, 
that the components near km s _1 toward £ Ori (where we have high S/N ECH-A spectra) also 
have relatively high excited state populations, and pressures significantly higher than those in the 
higher column density components near +25 km s _1 along that line of sight (Welty et al. 1999c). 
Higher resolution spectra of N I, O I, Na I, and Mg I A2853 would also be useful for more precise 
characterization of the individual WLV components. 



4. Predominantly Ionized Gas 

4.1. Ionized Gas Component Structure 

In addition to the strong absorption at low velocities, the strongest lines in the UV spectra 
also show evidence for the presence of additional lower column density components not seen in the 
optical lines — at low, intermediate, and high velocities (Figures § and ||). As discussed below, 
the gas at high and intermediate velocities is predominantly ionized; as noted above, there is 



The column densities of C II* for several of the WLV components are rather uncertain, due to blending with the 
weaker member of the A1335 doublet for the SLV components 17S and 18S, whose C II* column densities may be 
significantly larger (as discussed above). 
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some ionized gas at low velocities in addition to the dominant neutral gas. The HV gas ("Orion's 
Cloak"), at -108 km s" 1 <; v & -83 km s -1 , is seen in the strong lines of Mg II, Si II, Si III, Al II, 
N II, and C II; weak absorption is also detected for C II*, Al III, and perhaps S II (Figure ||). Four 
components (1H-4H) are present in the echelle profiles of C II, Si II, and Mg II. The IV gas, at 
—43 km s" 1 <^ v <^ —4 km s _1 (components 51-91), is seen most prominently in the strong lines 



of Si III, C II, Mg II, and (probably) N II. The component parameters listed in Table 10 for the 
generally weaker IV and HV components for several species were derived via a simultaneous fit to 
the echelle profiles of the stronger lines. That same component structure was used to fit the weak 
HV and/or IV absorption from Al III, S II, and C II*. For the high and intermediate velocity gas 
observed with G160M, we adopted the component structure determined from the echelle profiles 



of C II, Mg II, and Si II. Table 11 gives the resulting total column densities for the ionized HV and 
IV gas. 

Unfortunately, however, we have no high resolution spectra to define the component structure 
for the ionized gas at low velocities (apart from the ECH-B spectra of Al III, which shows weak 
components at —6 and +20 km s _1 ). This lack of component information particularly hinders our 
analysis of the strong N II A1083 and Si III A1206 lines. For Si III, we may obtain a lower limit to 
the IV + LV column density, log[A r (Si III)] > 13.4 cm" 2 , by integrating the "apparent" optical 
depth over the line profile (Hobbs 1971; Savage & Sembach 1991), but similar limits for N II are 
uncertain, because the true background level in the Copernicus spectrum is poorly known. We 
attempted to fit the Si III and N II absorption by using the HV and IV component structure found 
for C II (which, like Si III and N II, seems to be the dominant ion in the ionized gas toward 23 
Ori) together with either a "minimal" low velocity component structure based on the components 
found for Al III and Si II* or the WLV + SLV structure found for C II. The resulting HV column 
densities for Si III and N II are fairly well determined; the IV column densities are somewhat 
less so. For both species, iV(WLV) appears to be greater than iV(SLV), but the total A r (WLV + 
SLV) is very uncertain, so we have listed only (fairly conservative) upper limits in Table |ll|. The 
relatively weak absorption at low velocities (—5 km s _1 <^ v <^ +6 km s _1 ) seen in the G160M 
spectra of Si II* and S III (blended with C I and Si II) and of the higher ions C IV, Si IV, and 
N V (Figure ||) can be acceptably fitted in each case with a single broad component (Table |l2]) 
- though differences in central velocity and the relatively large 6-values suggest that additional 
substructure may be present. The weak absorption from the higher ions was not detected in the 
Copernicus spectra reported by Cowie et al. (1979), due primarily to the lower S/Ns achieved for 
those earlier spectra. In principle, P III could be present in the wing of the strong C II A1334 
line, but uncertainties as to the appropriate velocity and line width [and uncertainties in the SLV 
N(C II)] make it difficult to determine reliable and interesting limits. 



4.2. High- Velocity Gas 

4-2.1. HV Gas Abundances and Depletions 
No absorption from O I or N I is seen at high velocities, so the HV gas is primarily ionized 



(Trapero et al. 1996). The limit on iV(H I) in Table 11 assumes that O I is dominant in any 
neutral gas [which appears to be valid whether or not the gas is in equilibrium (e.g., Schmutzler 
& Tscharnuter 1993)] and that O is undepleted there. We do not detect high-velocity absorption 
from the higher ions Si IV and C IV, suggesting that only a narrow range of ionization states is 



-22- 



significantly populated for any element. Even so, determination of depletions is complicated by 
our inability to observe all the potentially significant ionization states of many of the elements. 
In principle, we would like to determine N{H II) by summing the contributions from the relevant 
ionization states of typically undepleted elements (e.g., C, N, and S), and then compare with 
the abundances of typically depleted elements (e.g., Si, Mg, Fe, and Al) to determine the level 
of depletion. We do not have data for C III, N III, Fe III, Mg III, or Al IV, however. Only for 
Si do we have essentially complete coverage — i.e., detections of Si II and Si III and significant 
limits on Si I and Si IV. The column densities of C II, N II, and Si II + Si III all yield estimates for 
log[iV(H II)] of about 17.7 cm -2 , assuming no depletion for C, N, and Si. While a similar estimate 
from N(S II) + N(S III) is consistent with that value, the uncertainties would permit a somewhat 
higher A r (H II), by perhaps a factor ~ 1.5, or 0.2 dex (assuming S undepleted)Q. Some C III and 
N III could thus be present, as found for the HV gas in several other lines of sight in Orion (Cowie 
et al. 1979; Welty et al. 1999c). In the HV gas toward C Ori, for example, Welty et al. (1999c) 
find N(C m)/N(C II) ~ 2.4 and JV(N m)/JV(N II) ~ 1.3. The ratios of iV(Al m)/JV(Al II), 
N(Si III) /iV(Si II), and N(S m)/JV(S II) are all significantly smaller in the HV gas toward 23 Ori 
than those found toward £ Ori, however, so that there is likely to be relatively less C III and N III 
toward 23 Ori. If log[iV(H II)] is greater than 17.7 cm -2 , then Si would be slightly depleted (but 
by less than 0.2 dex); likewise, C could not be depleted by more than 0.2 dex [and N(C III) must 
be less than N(C II)]. The ratio [iV(Al II) + N(A\ IIl)]/N(C II), however, is less than half the solar 
value, so that Al appears to be slightly depleted even for log[iV(H II)] ~ 17.7 cm" 2 . (The relatively 
low ratios of iV(Si III)/A r (Si II) and iV(Al IIl)/iV(Al II) suggest that there is unlikely to be an 
appreciable amount of Al IV present.) The depletion of Fe cannot be reliably estimated, since 
comparisons with the HV gas toward £ Ori suggest that a substantial fraction of the gas phase Fe 
may be in Fe III (Welty et al. 1999c). Slight differences in various ion ratios for the individual 
HV components may be due to differences in ionization for those individual components, as C II 
appears to be the dominant ion for C, but Si II and Mg II are not the dominant ions of Si and Mg. 

4-2.2. HV Gas Physical Conditions 

Cowie et al. (1979) ascribed the high velocity gas seen with Copernicus toward a number of 
Orion stars (including 23 Ori) to a radiative shock due to one or more recent supernovae in that 
region. They estimated a radius of about 120 pc (i.e., outside the expanding neutral shell) and 
an age of about 3 x 10 5 yr for this apparently large scale feature, as well as a pre-shock ambient 
density of about 3 x 10~ 3 cm~ 3 . Higher resolution GHRS spectra covering additional species have 
enabled a more detailed understanding of the composition and properties of this high-velocity gas 
(Trapero et al. 1996; Welty et al. 1999c). Toward 23 Ori, the newly acquired ECH-A spectra of 
C II, C II*, and Si II have yielded estimates for the temperatures and electron densities in the 
individual HV components. 

If the HV gas were in collisional ionization equilibrium, the observed ratios of various ions 
would be roughly consistent with an equilibrium temperature of about 25,000 K (using, for 
example, the calculations of Sutherland & Dopita 1993). The widths of the individual components 



10 We have combined our G160M spectra of the S III A1190 line with similar data from the HST archive to obtain 
a smaller upper limit on N(S III) in the HV gas than the one given by Trapero et al. (1996). 
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seen in Mg II, Si II, and C II, however, yield firm upper limits T max <J 14,000 K in all cases, and 
comparisons of 6(C II) vs. b(Mg II) and 6(Si II) suggest that the temperatures of the individual 
components are in fact 6,000-12,000 K (Table |l3|). The gas is therefore not in equilibrium, and 
has cooled more rapidly than it could recombine, presumably after having been shocked. We 
have found, however, that shock models are not entirely successful in reproducing the column 
densities, ion ratios, and line profiles observed for the HV gas in the Orion region (Welty et al. 
1999c). For example, models with u s hock ~ 100 km s _1 generally produce significantly more HV 
Si IV than the limits observed toward both 23 Ori and £ Ori. Furthermore, models of cooling gas 
(e.g., Schmutzler & Tscharnuter 1993) appear to predict significantly lower ionization for T < 10 
K than is observed for the HV gas in those two sightlines — though including the effects of an 
external radiation field may yield better agreement. 

In this ionized HV gas, the excited state C II* is populated primarily via collisions with 
electrons, and the ratio iV(C II*)/iV(C II) ~ 0.03 (for the three strongest components) can be used 
to estimate n e (and thus «h) from the equation of excitation equilibrium, since the temperature 
is known (T ~ 8000±2000 K): n e = 0.186T a5 [iV(C n*)/JV(C II)]. The n e = n H derived for the 
three main HV components lies in the range 0.4-0.5 cm -3 , about 150 times the pre-shock density 
estimated by Cowie et al. (1979). The thermal pressure, as measured by log(2n e T), is thus 3.7-4.0 
cm -3 K in each component. Estimating iV(H II) from iV(C II), assuming no depletion of C and 
no appreciable C III [neither would increase N(H II) by more than a factor 1.5], the individual 
component thicknesses range from about 0.005 pc to about 0.12 pc. 

4.3. Intermediate- Velocity Gas 

4-3.1. IV Gas Abundances and Depletions 

As for the HV gas, O I and N I are not detected at intermediate velocities (the G160M 
profiles of the strong A1302 and A1200 lines can be completely accounted for by the WLV and SLV 
components), and iV(Si III) > iV(Si II), indicating that the IV gas also is predominantly ionized. 
Some of the absorption from the higher ions S III, Si IV, C IV, and N V may also be associated 
with the IV gas, though lower ionization states are dominant (Table |ll]). From N(Si II) + 
iV(Si III) and iV(S II), we estimate that log[iV(H II)] ^ 17.7 cm -2 for the IV components; the less 
certain iV(N II) suggests that log[iV(H II)] could be as high as about 18.0 cm~ 2 . Since the iV(H II) 
estimated from C II alone is smaller by 0.2-0.5 dex than the above estimates, there may be some 
C III in the IV gas, as found toward several other Orion stars (Cowie et al. 1979). Comparisons of 
various ion ratios among those Orion lines of sight suggests, however, that iV(C III) is probably 
less than N{C II) in the IV gas toward 23 Ori — so that iV(N II) may be overestimated in 
Table |ll|. The ratio [N(A\ II)+iV(Al IIl)]/[iV(Si Il)+A r (Si III)] is similar to the value found for the 
HV gas, suggesting that Al is comparably depleted in the IV gas [unless there is an appreciable 
amount of Al in (unobservable) Al IV]. Higher resolution spectra of N II, Si III, and S III, as well 
as new spectra of C III, N III, and Fe III, would enable more reliable estimates for the abundances 
and depletions in the IV gas. 
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4-3.2. IV Gas Physical Conditions 

Ionized gas at intermediate velocities is seen toward a number of stars in the Orion region 
(Cowie et al. 1979), suggesting that the IV gas is a large-scale, regional feature — not intimately 
associated with the individual stars. Cowie et al. found the IV gas toward several of the Orion 
stars to be somewhat more highly ionized than the HV gas and so concluded that the IV gas is 
photoionized by (but at some distance from) the Orion association stars. In the IV gas toward 23 
Ori, both the slightly higher N(A\ III)/iV(Al II) and N(Si III)/iV(Si II) ratios (relative to those in 



the HV clouds; Table 11) and the (possible) associated absorption from C III, S III, Si IV, C IV, 
and N V would be consistent with somewhat higher ionization. The 6-values adopted for most of 
the IV C II, Mg II, and Si II components suggest, however, that the temperatures are comparable 
to those in the HV gas. Analysis of the C II fine structure excitation, assuming T ~ 8000 K, yields 
n e ~ 1.5-5.0 cm -3 , and thus log(2n e T) ~ 4.4-4.9 cm _3 K, for several of the IV components - 
somewhat higher densities and thermal pressures than the values found for the HV gas (Table |i~3| ). 
Several of the IV components have significantly lower iV(H II), however, and so are much thinner 

0.001 pc) than the HV components — suggesting that they may arise in thin filaments or 
sheets. 



4.4. Ionized Gas at Low Velocity 

The strong absorption seen for N II and Si III and the weak, broad absorption due to S III, 
C IV, Si IV, and N V indicate that there is some ionized gas at low velocities (Figure ^ and 
Tables |ll| and 12). The somewhat narrower Si II* absorption (centered at about +12 km s _1 ) and 



the low-velocity C II* absorption may be largely due to neutral gas, however (§§ 3.3.2 and 



Because some of the broad, unresolved absorption from the higher ions may be associated with 



the IV gas, we have listed the total high ion column densities in a separate column in Table |11. 
As discussed above, the low-velocity column densities for N II and Si III are rather uncertain, but 
the limit for N II suggests that log[iV(H II)] <^ 19.3 cm~ 2 . The observed column density of S III 
provides a lower limit of log[iV(H II)] ;> 18.2 cm~ 2 (though some of this could be at intermediate 
velocities). The 6-values found from single component fits to the high ion profiles — 16-17 km s _1 
for S III and Si IV, 23-24 km s" 1 for C IV and N V) — provide upper limits for the temperature of 
4-5 x 10 5 K. The b- value for the +20 km s _1 component of Al III, however, implies T <^ 3 x 10 4 
K, and it seems likely that the low-velocity ionized gas is comprised of several components. There 
is no obvious correlation between central velocity and ionization potential Xion f° r these high ions. 
We will briefly examine several possible origins for the high ions in Section pT 



5. Discussion 

5.1. The Orion Region 

General surveys of interstellar Na I and Ca II absorption have revealed complex multi- 
component absorption-line profiles toward a number of the brighter stars in the Orion-Eridanus 
region (Hobbs 1969, 1978b; Marschall & Hobbs 1973; Beintema 1975; Welty et al. 1994, 1996). 
Frisch et al. (1990) obtained Ca II and/or Na I spectra of a number of stars in that region in 
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order to locate the various velocity components along the lines of sight. We will provide a brief 
update of the latter work, making use of additional (as yet unpublished) optical spectra and the 
parallaxes from Hipparcos. In the southern part of this region (within about 5° of A Eri), weak 
Ca II absorption at +20 km s _1 ^ v ^ +30 km s _1 is seen in the spectra of several stars at 
distances d ~ 170 pc, but is not detected (at the level of several mA) toward stars at d ^ 100 
pc. Nearer to 23 Ori, several weak components are seen in both Na I and Ca II toward 7 Ori (2?9 
N of 23 Ori), at d ~ 75 pc, but no Na I is detected toward HD 35134 (0?8 S), at d ~ 130 pc. 
Stronger absorption at those velocities, with significantly larger iV(Na l)/iV(Ca II) ratios, is seen 
over most of the region toward stars beyond about 270 pc — including HD 34959, 25 Ori, 30 Ori, 
and 33 Ori, all within 2° of 23 Ori. Weaker absorption at —10 km s _1 <J v <J +20 km s^ 1 , with 
lower iV(Na I)/iV(Ca II) ratios, is also seen toward many stars beyond that point. While there are 
significant differences in the column densities and velocities of individual components seen toward 
the stars beyond 270 pc, the overall general similarities among the profiles suggest that the bulk 
of the observed gas is not closely associated with the individual stars, but is likely to be located in 
the foreground, at 130 pc <^ d <^ 270 pc. The observed differences are not unexpected, in light of 
the variations commonly seen on much smaller spatial scales (Watson & Meyer 1996) and of the 
energetic activity in the region. 

Comparisons among the 21 cm, Ha, soft X-ray, and 100 fim emission maps of the Orion- 
Eridanus region have suggested the presence of a large, expanding shell, predominantly neutral 
but with an inner ionized zone, encompassing a bubble of hot gas at log(T) ~ 6.2 K (Heiles 
1976; Reynolds & Ogden 1979; Burrows et al 1993; Brown et al. 1995; Snowden et al. 1995; see 
also Cowie et al. 1979). If that picture is accurate, 23 Ori would lie within the bubble, near its 
boundary at low Galactic latitudes. The SLV (and perhaps the WLV) components could arise in 
the expanding shell, based on the distance estimates given above and on the general agreement 
in velocity with the shell's 21 cm emission near the line of sight to 23 Ori (Brown et al. 1995). 
The WLV gas may represent the more "active" part of the shell, in view of its higher temperature 
and thermal pressure, more negative velocities, and less severe depletions. The total H I column 
density, 5.5 x 10 20 cm 2 , is also very similar to that inferred for the shell (Snowden et al. 1995). 
The ionized low- velocity gas could then represent the inner ionized region of the shell, estimated 
to contain about 10-15% of the total hydrogen in the shell (Reynolds & Ogden 1979). If 23 Ori 
itself is immersed in the hot, low density gas contained within the shell, we would not expect a 
localized, associated circumstellar H II region — consistent with the lack of absorption due to 
ionized gas observed within 10-15 km s _1 of the stellar radial velocity. (In any case, a Bl V star 
would typically have only a very small associated H II region.) 

If the gas seen in absorption toward 23 Ori (and toward other Orion OBI stars) is not, in 
general, closely associated with the individual background stars, then the radiation field incident 
on the clouds will not be dominated by the flux from those individual stars, but will be due to 
the ensemble of stars in the region. Wall et al. (1996) have estimated the radiation field in the 
Orion region by assuming all the stars in the Orion OBI and A Ori associations to be located at a 
distance of 450+10 pc, with some allowance for extinction effects. In that model, the field along 
the line of sight to 23 Ori is roughly twice the "average" interstellar field of Mathis, Mezger, & 
Panagia (1983) at d ~ 450 pc — and presumably less than that for smaller distances. Another 
estimate for the radiation field may be obtained by considering the relationship between N(H.2) 
and iV(H I), as a function of temperature, density, and radiation field (Reach, Koo, & Heiles 
1994). Using equation (3) of Reach et al., with the N(H 2 ) = 2 x 10 18 cnT 2 and N(H I) = 5.1 x 



-26- 



10 20 cm -2 observed for the SLV clouds and the T = 100 K and bh ~ 10-15 cm -3 derived above 
(§ |3.3| ), we would find the radiation field to be enhanced by a factor \ ~ 2.5-3.8 over that average 
interstellar field. That equation, however, is applicable only if the H2 has reached its equilibrium 
abundance. Consideration of the H2 formation time (Reach et al. 1994), for the T and tt-h 
obtained for the SLV clouds, yields if or m ~ 7 x 10 8 yr (i.e., much longer than the apparent age of 
the Orion-Eridanus bubble) — suggesting that the H2 has not reached equilibrium. If the eventual 
equilibrium iV(H2) is higher than the currently observed value, then the true radiation field is 
likely to be somewhat weaker, since 7V(H2) oc x~ 2 - We have assumed in discussions above that 
the radiation field incident on the clouds toward 23 Ori can be approximated by the "average" 
WJ1 interstellar field defined by de Boer, Koppenaal, & Pottasch (1973), which agrees with the 
Mathis et al. field to within 20% below 1300 A. In light of the above estimates, the radiation field 
could be somewhat stronger, but probably not by more than a factor 2. The field could also be 
somewhat "harder" , in view of the number of bright O stars in the region and the shallow far-UV 
extinction observed toward some Orion stars (e.g., Fitzpatrick & Massa 1990). Toward Orion, the 
field is unlikely to be weaker than the WJ1 field (except in significantly shielded regions). 



5.2. Ionization and Electron Densities in the SLV Gas 

Because the SLV components exhibit detectable absorption from an unusually large number 
of trace neutral species (including Si I, Ca I, Fe I, and perhaps Mn I), we have many independent 
diagnostics for the electron density in the SLV gas. As described above (§ |3.3|) , the estimates for n e 
resulting from the assumption of simple photoionization and radiative recombination are generally 
relatively large (implying n e /nn ~ 0.01), but cover a wide range, from about 0.04 cm -3 to about 
0.95 cm -3 (Table |8|). Detailed studies of other lines of sight show analogous, seemingly systematic 
offsets among the n e inferred from different diagnostics — but with subtle variations among 
the element-to-element offsets for the different lines of sight (Cardelli, Sembach, & Savage 1995; 
Fitzpatrick 1996, 1998; Snow et al. 1996; Fitzpatrick & Spitzer 1997; Welty et al. 1999a, 1999b). 
Several of those studies also find evidence for some ionization of H in predominantly warm, diffuse, 
neutral gas. The relatively high n e /nn in the SLV components toward 23 Ori, however, apparently 
occurs in colder, denser gas. Accurate, consistent estimates for both n e and n e /nn are necessary 
to understand what processes are responsible for the ionization (simple photoionization, cosmic 
ray ionization, charge exchange processes, decaying neutrinos, ...) and recombination (radiative, 
dielectronic, charge exchange, ...) in these interstellar clouds. In studies of diffuse clouds, the 
fractional ionization is often assumed to be ~ 3 x 10~ 4 (due to photoionization of C and other 
heavy elements) in order to estimate bh from n e — but that assumption may be invalid in many 
cases. We will discuss these issues in more depth in a separate paper (Welty et al. 1999b), but 
briefly present some of the possibly relevant processes here as they apply to the SLV gas toward 
23 Ori. 

• The neutral species Ca I and Mg I can be enhanced via dielectronic recombination, for 
temperatures ;> 3000 K and ;> 6000 K, respectively — which would then yield apparent 
overestimates of n e assuming only radiative recombination. As discussed above (§ |3.3|) , 
however, the temperatures in the SLV components are likely to be ~ 100 K. While the 
iV(Ca l)/iV(Ca II) ratio often yields much higher estimates for n e than do other neutral-first 
ion ratios, the Ca I seems to arise predominantly in cool gas in the Galactic ISM, so its 
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apparent enhancement does not, in general, seem to be due to dielectronic recombination 
(Welty et al. 1999b). 

• The derived column densities could be in error, due either to errors in the /-values or to 
problems with the profile-fitting analysis. In many cases, however, the column densities were 
derived from several lines of different strength — including some very weak lines — which 
yielded consistent results. In most cases, the nominal uncertainties in the /-values are less 
than about 0.1 dex (Morton 1991), and a number of the /-values have more recent, more 
accurate determinations (see the Appendix). 

• The photoionization and/or recombination rates for some species could be in error, due to 
inaccurate cross-sections or to differences between the true and adopted radiation fields. 
A recent estimate of the rates for the Fe I— Fe II equilibrium, for example, suggests that 
T/a is roughly a factor 2 higher than previously thought — apparently due largely to a 
more accurate and consistent accounting for the strong resonances in the photoionization 
cross-section of Fe I (Nahar, Bautista, & Pradhan 1997; see also Bautista, Romano, & 
Pradhan 1998). A higher T/a would increase the (currently low) n e inferred from the 
N(Fe l)/A r (Fe II) ratio. For most of the species re-considered by Bautista et al., the new 
T's are larger than previously estimated, but new, consistent a's are needed in order to 
obtain revised T/a ratios. New estimates for a have been determined for C and Si (Nahar 
& Pradhan 1997), yielding T/a ratios slightly smaller than previous values, as both T and 
a increased by similar amounts. Self-consistent re- appraisals of the rates for other elements 
would be very useful — though the net effect on the T/a ratios may generally be minor, if the 
issue is primarily one of accounting for the resonances in the photoionization cross-section. 

In H I regions, the photoionization rates for species with relatively large ionization potentials 
(C I, S I, P I, CI I, Ca II) depend on the (poorly known) radiation field in the narrow interval 
between xion and the Lyman limit at 13.6 eV. The rates for those species would thus be 
somewhat higher if the field in the Orion region is harder than the WJ1 field — as illustrated 
in Table || for the Draine (1978) radiation field, which is stronger than the WJ1 field below 
about 1400A. The rates for other neutrals with lower Xion ar e dominated by the field at 
longer wavelengths, and thus depend much less on the choice of radiation field. The neutral 
species C I and S I do appear to be relatively enhanced for lines of sight characterized by 
moderate reddening [E(B-V) <; 0.3] and steep far-UV extinction, where the far-UV field is 
preferentially suppressed (Welty 1989; Welty et al. 1991; see also Jenkins &: Shaya 1979). 
The expected magnitude of this effect depends on the grain scattering properties, and will 
be smaller if the grains are very forward scattering in the far UV. This effect should be 
insignificant for the clouds toward 23 Ori, however, since the color excess is small and the 
far-UV extinction does not appear to be steeper than "average" (Savage et al. 1985; Papaj, 
Wegner, & Krelowski 1991). Furthermore, differences in the radiation field in the far-UV 
would tend to affect C and S similarly (relative to other elements) — which does not appear 
to be the case toward 23 Ori — though other processes may affect both (see below). 

• In principle, a dependence of depletion on local density nn could produce differences in 
inferred n e between mildly and severely depleted elements. If the density of trace neutral 
species is proportional to (assuming, for example, n e oc hh; though see Welty et al. 
1999b), then the trace neutrals will be more concentrated in the denser regions of a cloud. 
(Note that this implies that all the n e estimates from trace/dominant ratios underestimate 
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the n e in the denser regions.) Integrated over the whole cloud, the ratios iV(X I)/iV(X II) 
for depleted elements will be dominated by the less dense regions, while the ratios for less 
depleted elements will be more heavily weighted by the denser regions — which could lead 
to systematic differences in derived n e , if n e depends (in some way) on ran- In the SLV 
clouds, however, there does not appear to be any relationship between elemental depletion 
Sx (as determined from the dominant species) and the n e inferred from iV(X l)/iV(X II). For 
example, the lowest n e are derived from S, Mn, and Fe — which have significantly different 
depletions [column 5 (n e ) versus column 6 ([X/C]n = Dx — Dc) in Table |j. 

• Fitzpatrick & Spitzer (1997) derived systematically higher n e from the iV(C l)/iV(C II) ratio 
(versus several other ratios) in the halo clouds toward HD 215733, under the assumption of 
simple photoionization and radiative recombination^. They conjectured that dissociative 
recombination of CH + might have enhanced the C I, though CH + has not been observed 
toward HD 215733. At T ~ 100 K, the radiative recombination rate for C II to C I is a 

~ 9 x 10~ 12 cm 3 s _1 (Nahar & Pradhan 1997) and the dissociative recombination rate for 
CH+ is (3 ~ 3.3 x 10~ 7 cmV 1 (Mitchell & McGowan 1978). For the SLV clouds toward 23 
Ori, log[iV(C II)] ~ 16.95 cm' 2 and log[iV(CH + )] ~ 13.06 cm" 2 . We would thus find that 
dissociative recombination should dominate radiative recombination by a factor of about 4.8 

- which would imply that the n e derived from C should be reduced to about 0.04 cm -3 , 
similar to the values derived from S, Mn, and Fe. That comparison assumes that the CH + 
coexists in the cool SLV clouds with the various atomic species. The CH + molecular ion, 
however, is rapidly destroyed (by H2 and by H atoms as well as by electrons) in cold, dense 
regions. Furthermore, if we fit the CH + profile (either ours or the higher resolution spectrum 
obtained by Crane et al. 1995) with two components corresponding to the two strongest 
SLV components, we obtain 6-values that are larger than those that would result from T ~ 
100 K and the turbulence implied by the SLV K I or Ca II line widths. It is thus likely that 
the CH + is not completely coextensive with the various atomic species in these relatively 
cold clouds, even though the component velocities are similar. The effect of dissociative 
recombination on the carbon ionization balance is thus somewhat uncertain. 

Given the current limit log[iV(SH + )] < 12.7 cm' 2 toward 23 Ori (Ma gnani & Salzer 1991), 
dissociative recombination could also enhance iV(S I), if the rate is of order 10 -7 cm 3 s -1 
(Fitzpatrick 1998). The S I column density is low, relative to the other trace neutral species, 
however, which may imply that the SH + abundance is significantly smaller than the current 
upper limit. 

• We have not included chlorine in Table |8|, because its ionization balance is dominated by 
a rapid reaction between CI II and H2, with rate constant k = 7 x 10~ 10 cm 3 s _1 , which 
ultimately leads to enhanced CI I (Jura 1974; Fehsenfeld & Ferguson 1974; Jura &: York 
1978). For the SLV clouds, if we take n(H 2 ) ~ iV(H 2 ) / 14 pc (likely a lower limit to the 
H-2 density), the CI II radiative recombination rate a = 8 x 10~ 12 cm 3 s _1 (Pequignot & 
Aldrovandi 1986), and n e ~ 0.15 cm" 3 , then [n(H 2 ) k] j [a n e ] ^ 30, so the reaction with H 2 
is much faster than radiative recombination. If we estimate the chlorine depletion as Dei ~ 



1 We note, however, that use of the revised /-value for the Mg II A1239 doublet (Fitzpatrick 1997; see the 
Appendix) brings the n e determined from C and Mg toward HD 215733 essentially into agreement. Furthermore, if 
Ca is more depleted than Fe in those clouds, as we find for the WLV clouds toward 23 Ori, then the n e determined 
from Ca II/Fe II may also be in agreement. 
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-0.4 dex (i.e., slightly less severe than the cold cloud value, as found for other lightly depleted 
elements in the SLV clouds), then log[A(Cl II)] ~ log[iV(Cl I)] ~ 13.3 cm" 2 , consistent with 
the very uncertain value given by Jenkins et al. (1986). The photoionization rate T for CI I 
can then be estimated from the equation T 7V(C1 1) = k n{B.2) N(C\ II). The resulting value 
is about twice the T = 1.6 x 1CP 11 s _1 given by Pequignot & Aldrovandi (1986) for the WJ1 
field — perhaps suggestive of a somewhat harder far-UV radiation field (xion = 12.97 eV for 
CI I). 

• Lepp et al. (1988) suggested that charge exchange between singly ionized species and either 
neutral or negatively charged large molecules (e.g., PAH's) might significantly increase 
the relative fractions of neutral species in interstellar clouds — thus leading to spuriously 
high n e inferred from simple ionization balance calculations. The computed neutral species 
enhancements, for models of the main clouds toward £ Per and £ Oph, range from factors 
of about 10 to 15 for Mg, Ca, and Fe, depending on the details of the cloud model (and on 
the rather uncertain properties of the PAH's), as the PAH abundance increases from 10~ n 
to 10~ 6 (relative to hydrogen). The enhancements appear to be about 30% smaller for C 
and S, presumably due to their larger ionization potentials; estimates were not given for 
other elements. Thus, while charge exchange with large molecules might provide a partial 
explanation for the apparently large overall n e (and n e /nn) in the SLV clouds, it would not 
appear to account for the large differences in n e inferred from different elements. 

• Pequignot & Aldrovandi (1986) suggested that charge exchange between neutral species and 
protons or other positive ions (X° + H+ - > X+ + H° + AE) could affect the abundances 
of the neutral species in cool interstellar clouds. They identified which reactions would be 
significant and made rough guesses as to the rates, based on the energy (AE) released by 
the different possible reactions. Where the fractional ionization is low [e.g., in the main ( 
Oph ("B") clouds considered by Pequignot & Aldrovandi], reactions with C + can dominate; 
where the fractional ionization is high (as appears to be the case in the SLV clouds toward 
23 Ori), reactions with protons should dominate, essentially to reduce some of the neutral 
species column densities. In the latter case, the species Al I, S I, Ti I, Mn I, Fe I, and Ni I 
should be most affected; for all these, Pequignot & Aldrovandi estimated a rate constant 

(3 ~ 3 x 10~ 9 cm 3 s -1 . [We note, however, that the relevant charge exchange rates have 
not been determined for most of the species in question for temperatures as low as 100 K 
(Kingdon k, Ferland 1996; Stancil et al. 1998).] Interestingly, S I, Mn I, and Fe I are the 
three apparently underabundant neutral species in the SLV clouds toward 23 Ori (and Ni I 
may also be underabundant). If we add a charge exchange term to the ionization balance 
equation, we have X°T + (3X°n p = an e X + (with the proton density n p = n e ). If we assume 
= 3 x 10~ 9 cm'V 1 and T = 100 K, then we infer n e = 0.06, 0.18, and 0.10 cm' 3 for S, 
Mn, and Fe, respectively — closer to the estimates derived from most other elements than 
the 0.04-0.05 cm~ 3 obtained with no charge exchange. If we instead fix n e = 0.15 cm" 3 (see 
below), then we would need rate coefficients f3 = 10.3, 2.8, and 3.5 x 10~ 9 cm 3 s _1 for S, 
Mn, and Fe to obtain that n e . Charge exchange between protons and neutral species thus 
appears to provide a possible explanation for the low n e derived from S, Mn, and Fe in the 
SLV clouds toward 23 Ori, if in fact n e /nn is as high as 0.01 (for protons to dominate) and 
if the rates are of order a few times 10~ 9 cm 3 s~ 1 . Calculations by Kimura et al. (1997) 
indicate, however, that the cross-section for charge transfer ionization of S I is smaller than 
the corresponding cross-section for Si I (Kimura et al. 1996) — suggesting that the rate 
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coefficient for S I is of order 10 11 cm 3 s 1 at 100 K, much smaller than would be needed 
to significantly affect the sulfur ionization balance (above; Pequignot & Aldrovandi 1986; 
Fitzpatrick & Spitzer 1997). 

The large range in n e inferred from different trace/dominant ratios in Table || and the 
uncertainties associated with the various additional processes that might affect the ionization 
balance for different elements imply that any estimate of n e will be somewhat uncertain. We have 
adopted n e ~ 0.15±0.05 cm -3 (from Na, Mg, Si, and K) as an average value for the SLV clouds 
toward 23 Ori [and recall that n e (17S) ~ 2.6 n e (18S)]. Charge exchange between neutral species 
and protons may account for the lower n e inferred here from S, Mn, and Fe under the assumption 
of simple photoionization and radiative recombination, though more accurate values for the 
charge exchange rates at low temperatures are needed to confirm this. A harder far-UV radiation 
field would lead to lower n e estimated from C, P, and S (assuming the WJ1 field); dissociative 
recombination of CH + , however, would increase the estimate from C. At the low temperatures 
found for the SLV clouds, dielectronic recombination should not affect the n e estimates from Mg 
and Ca. Errors in the adopted photoionization and/or recombination rates (e.g., for Fe) may also 
contribute to the range in n e derived from different diagnostic ratios. The true n e could be higher 
if the actual radiation field is higher than the adopted "average" Galactic field; the actual field is 
unlikely to be lower than average. The true n e could be lower, by perhaps an order of magnitude, 
however, if charge exchange between singly ionized species and large molecules is important at 
the moderate densities found for the SLV clouds. A smaller n e would be more consistent with 
estimates based on the iV(C II*)/iV(C II) ratio, as discussed above (§ |3.3.3| ; see also Fitzpatrick 
1998). Even if that is the case, the n e /nn ratio would still imply some partial ionization of 
hydrogen in the SLV clouds. 

For other lines of sight, we would emphasize that inferred n e are likely to be very unreliable, 
apart from a detailed understanding of the physical conditions and processes affecting the 
ionization balance. In many cases, for example, the n e obtained from iV(Ca l)/A(Ca II) are 
higher than the n e inferred from other such ratios — by an order of magnitude or more. Such a 
systematic effect could be quite significant, as the iV(Ca l)/iV(Ca II) ratio has been extensively 
used both to calculate electron densities and to infer abundances for elements (e.g., Na, K, Li) 
whose dominant ionization states cannot be observed. Estimates of relative depletions based on 
ratios of trace neutral species (e.g., Snow 1984; Snow et al. 1996) will also be affected by any 
additional processes selectively affecting the ionization balance. 

There are several possible explanations for the partial ionization of the hydrogen in the 
predominantly neutral SLV gas. The enhanced soft X-ray flux observed from the Orion-Eridanus 
region (Burrows et al. 1993), which may be due to thermal emission from hot gas [log(T) ~ 6.2 
K] produced by Orion OBI association stellar winds and/or supernovae, could ionize some of the 
H I in the surrounding neutral clouds. Snowden et al. (1995) estimated that an intervening H I 
column density comparable to that observed toward 23 Ori could account for the smaller ratio of 
b/4 keV to 3 /4 keV flux toward Orion, relative to the ratio seen in the part of the bubble farther 
from the Galactic plane in Eridanus. We may use equation 5-32 in Spitzer (1978) to estimate the 
ionization rate required to produce the observed fractional ionization in the SLV gas: 
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Here a^ 2 \ the hydrogen recombination rate (excluding recombinations to the lowest level), is ~ 
7 x 10~ 12 s _1 for T ~ 100 K, and we set nj, the density of positive ions heavier than He, to 2 
x 10~ 3 cm~ 3 , based on the observed gas-phase heavy element abundances (primarily C) and the 
derived n(H I) ~ 10 cm~ 3 . To obtain the inferred fractional ionization n e /rin ~ 0.01, we thus need 
an ionization rate Ch ~ 7 x 10 -15 s _1 . While this rate is two orders of magnitude larger than the 
cosmic ray ionization rate usually assumed for the local ISM, it is only a factor of 2 larger than 
the rate inferred for some of the clouds observed toward HD 93521 and HD 149881 (Spitzer & 
Fitzpatrick 1993, 1995), and some enhancement might be expected in the Orion region, in view of 
the X-ray observations noted above. If charge exchange with large molecules has enhanced the 
neutral species by a factor of 10, however, the inferred ionization fraction would be 0.001, requiring 
an ionization rate lower by a factor of 100, since Ch ~ (a^ n 2 ) / tih- Quantitative estimates of 
the X-ray flux — and the predicted ionization of hydrogen — would be valuable. 

If the SLV components arise in the expanding shell around the Orion-Eridanus bubble, the 
large apparent fractional ionization might also be due to partial mixing of the dominant neutral 
gas with the adjacent (interior) layer of ionized gas, which Reynolds & Ogden (1979) estimate 
comprises about 10% of the material in the shell. As noted above, some ionized gas is present 
at similar velocities, and such a mixture would have only minor effects on the abundances and 
depletions in the SLV gas. This possibility is intriguing, in light of the relatively high abundances 
of CH and CH + , as turbulent boundary layers have been proposed as one possible source of CH + 
(e.g., Duley et al. 1992, and references therein). Duley et al. present diagnostic ratios (e.g., 
ratios of the densities in several H2 rotational levels to that of CH + ) for steady-state models 
characterized by several different temperatures, densities, and radiation fields. The models with T 
= 2000 K, n = 20 cm" 3 , ra(CH+) ~ 1.5 x 10~ 4 cm , and radiation field 1-3 times the average IS 
field seem to agree best with the ratios of CH and of the J=3-5 levels of H2, with respect to CH + , 
toward 23 Ori. If those models are appropriate for the 23 Ori line of sight, the CH, CH + , and the 
higher rotational levels of H2 might thus arise in a thin (~ 0.02 pc), warm layer adjacent to the 
colder SLV gas. Observations of OH absorption toward 23 Ori could provide a stringent test for 
the Duley et al. models. 

Either of these possible explanations for increased n e /nn toward 23 Ori would also apply 
to the corresponding strongest neutral components toward other stars in the Orion region. 
Preliminary analyses of the strongest components at low velocity toward C Ori seem to indicate 
that both n e and the fractional ionization are somewhat lower than in the SLV clouds toward 23 
Ori, but that the pattern of neutral/first ion ratios for different elements (and thus the inferred 
relative n e ) appears to be very similar. 

5.3. Elemental Abundances and Depletions 

5.3.1. Previous Studies 

Spectra obtained with Copernicus and IUE provided total line-of-sight abundances and 
depletions for a number of elements (e.g., Jenkins 1987; de Boer et al. 1987), but the relatively low 
resolution of those spectra has generally precluded studies of individual interstellar clouds. The 
average line-of-sight depletions typically exhibit a fairly smooth variation with average line-of-sight 
hydrogen density < «h > = N(R)/d (though with some scatter), becoming increasingly severe 
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for larger values of < bh >. Similar trends have been noted versus the fraction of hydrogen in 
molecular form /(H2), which may be a better indicator of local physical conditions than < bh > 
(e.g., Cardelli 1994). Spitzer (1985) proposed a simple, idealized model, in which a given line of 
sight intercepts a mixture of several distinct cloud types, each having a characteristic overall level 
and pattern of depletions. In that model, warm, low density clouds have relatively mild depletions, 
while colder clouds (both "standard" and "large") are characterized by more severe depletions. 
Jenkins et al. (1986) used that model to interpret a large body of Copernicus abundance data 
for Mg, P, CI, Mn, Fe, Cu, and Ni [see Jenkins (1987) for additional elements], and derived 
representative depletion values for the cold and warm clouds. We have preserved some of this 
framework in the "cold cloud," "warm cloud," and "halo cloud" depletions listed in Table but 
we view the values as representing a likely continuum of depletions encountered in individual 
clouds characterized by different temperatures, local densities, and histories. Joseph (1988, 1993) 
took the "best" of the Jenkins et al. (1986) data, and found the depletions of Mg, P, CI, Mn, 
and Fe to be linearly correlated — contrary to expectations, based on ideas of grain processing 
in the ISM, that some elements should deplete more readily than others (e.g., Barlow 1978). 
Joseph conjectured that the relative depletions were fixed at grain formation, and that subsequent 
processing in the ISM (accretion, sputtering, etc.) preserved those relative depletions. Several 
attempts to detect preferential depletions, by comparing several lines of sight in a common region 
but with different overall levels of depletion (Snow & Jenkins 1980) or by using trace neutral 
species to probe denser regions (Snow 1984; Snow, Joseph, & Meyer 1986) likewise failed to reveal 
any strong, consistent signature of differential depletion. In retrospect, several factors — blending 
and possible saturation effects in the relatively low-resolution UV spectra, a somewhat limited 
range in overall depletions for the elements considered, and other processes affecting the neutral 
species — may have conspired to hide element-to-element differences in depletion behavior. We 
note that higher resolution optical spectra of Ti II (e.g., Stokes 1978) had already suggested that 
the Ti depletion exhibits a significantly larger range (compared to less severely depleted elements) 
for lines of sight largely overlapping those used in the UV surveys. 

The higher spectral resolution achievable with the HST GHRS has brought us closer to the 
goal of determining the abundances, depletions, and physical properties — and their possible 
interrelationships — for individual interstellar clouds. Savage &: Sembach (1996) have summarized 
the results of a number of GHRS studies of interstellar abundances and depletions; Fitzpatrick 
& Spitzer (1997) presented a detailed study of the individual components found toward the 
halo star HD 215733. In one of the first individual-line-of-sight abundance studies using spectra 
from GHRS, Savage et al. (1992) compared the abundances and depletions for two component 
groups (the aforementioned groups A and B) seen toward £ Oph, and found that the differences 
in depletion between the two groups were larger for the more severely depleted elements Fe, 
Cr, Ni, and Ti than for the less depleted O, Mg, Cu, and Mn (see Figure ^) — i.e., differential 
depletion. Sofia et al. (1994), incorporating data for several additional lines of sight, presented 
evidence for the preferential depletion of Mg and Fe relative to Si, developed a core-mantle model 
for the dust grains, and conjectured that the total (gas + dust) abundances in the ISM might 
be sub-solar. Snow et al. (1996) inferred preferential depletion of some elements in the denser 
regions of the clouds toward HD 154368 from comparisons of trace neutral and dominant singly 
ionized species, but ionization balance uncertainties (§ |5.2| ) may invalidate those conclusions. 
Many of the other lines of sight observed with the GHRS (including a number toward stars in 
the Galactic halo) sample primarily warm, fairly diffuse gas, and have provided constraints on 
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the minimum depletions in the diffuse ISM, due perhaps to a population of relatively hardy grain 
cores (Fitzpatrick 1996). 

5.3.2. 23 Ori vs. ( Oph — Density-Dependent Depletions? 

Comparisons among the depletions obtained for the 23 Ori SLV components, those found for 
the £ Oph B components, and the values adopted for "typical" cold, dense clouds suggest that a 
range of depletions may be present for relatively cold, neutral clouds and that the differences in 
depletion may be related to differences in local physical conditions. Except for the very lightly 
depleted C, N, O, S, Sn, and Kr, the SLV depletions are 0.3-0.6 dex less severe than the adopted 
cold cloud values (relative to H; i.e., [X/H]), whereas the ( Oph B depletions are comparable to or 
slightly more severe than the cold cloud values. This is the case for lightly depleted elements (P, 
Zn), for moderately depleted elements (Mg, Si, Mn, Cu, Ge, Ga), and for more severely depleted 
elements (Al, Ca, Ti, Cr, Fe, and Ni). We note, however, that since Zn also is less depleted than 
in the cold cloud pattern, the SLV depletions of other elements, relative to Zn (i.e., [X/Zn]), are 
more similar to the cold cloud pattern. Since the SLV depletions of C, N, O, and Kr are consistent 
with those obtained for other lines of sight [essentially independent of < wh > or /(H2)], and 
since several independent determinations find the total log[iV(H I)] less than or equal to the 
adopted 20.74 cm~ 2 , it is unlikely that we have significantly underestimated the amount of neutral 
H in the SLV components. Furthermore, there is apparently only a small amount of ionized H 
at those velocities. Finally, the individual component values of iV(Fe Il)/iV(Zn II), for the two 
strongest SLV components, differ by only about 0.2 dex. The differences in depletion between the 
23 Ori SLV components and the ( Oph B components thus do not appear to be due to errors in 
abundances, mixtures of very disparate components, or contamination from ionized gas — but 
instead appear to reflect a range of at least a factor 2 in the depletions among different cold clouds 
(cf. also the comparisons between the depletions toward £ Per and £ Oph in Savage et al. 1992 
and toward 1 Sco and £ Oph in Welty et al. 1995). 

The relationship between average line-of-sight depletions and < wh > recognized in the early 
UV data was suggestive of a dependence of depletions on local hydrogen density. Joseph et al. 
(1986) analyzed /[/i? spectra for several lines of sight showing strong CN absorption (presumably 
indicative of dense gas), and found some evidence for preferential depletion of Mn and perhaps 
Cr, but a number of the abundances are rather uncertain. Cardelli, Federman, & Smith (1991) 
compared Ca I abundances for lines of sight with different iV(CN)/iV(CH) ratios, and concluded 
that the Ca depletion is proportional to n^ 3 . All those studies, however, were based on integrated 
line-of-sight abundances, and actual values of «h were not known in most cases. It is therefore 
intriguing that the densities derived for the SLV components toward 23 Ori (nn ~ 10-15 cm~ 2 ) 
are at least an order of magnitude smaller than those found for the £ Oph B components, where 
the depletions are more severe by a factor of 2-4. For the main component (group) B toward 
C Oph, Lambert et al. (1994) obtained wh ~ 100-200 cm~ 2 from C I fine-structure excitation, 
Federman et al. (1994) found nn ~ 400 cm -2 from an analysis of the carbon chemistry, and Lepp 
et al. (1988) used central densities of 350 and 800 cm' 2 in their models. On the other hand, if the 
SLV components are located in a shell of gas swept up by the expanding Orion-Eridanus bubble, 
we might expect the depletions to be somewhat less severe — due either to partial erosion of the 
associated grains or to incorporation of lower density gas with milder depletions. Detailed studies 



-34- 



of other lines of sight are needed to determine whether there is any general relationship between 
the depletions and the local densities for individual cold clouds. 

Snow (1984; see also Snow et al. 1986, 1996) investigated the possible dependence of depletion 
on density by comparing the relative depletions inferred from ratios of trace, neutral species 
(which may be concentrated in the denser regions) with those obtained from the presumably more 
broadly distributed dominant species. Significant differences in the relative depletions inferred 
from the trace and dominant species were found only toward £ Oph and HD 154368. As noted 
above (§ |5.2|) , however, such comparisons assume that photoionization and radiative recombination 
dominate the ionization equilibria for the various elements considered. In the last two columns of 
Table ||, we list relative depletions for a number of elements — inferred both from the dominant 
first ions and from the trace neutral species. In both cases, the depletions are given relative to 
that of carbon, whose depletion appears to be relatively constant over wide ranges in both < ran > 
and /(H2) (Sofia et al. 1997). For all elements except Ca, the relative depletions inferred from the 
neutral species are more severe than those determined from the dominant species. The differences 
are less than a factor 2 for Na, Mg, and K, but may be as much as factors of 5-7 for S, Mn, and 
Fe — corresponding to the differences in n e inferred from those elements. On the other hand, if 
the depletions are instead computed relative to that of sulfur (which is typically undepleted), then 
in many cases the relative depletions inferred from the neutral species are less severe than those 
determined from the dominant species. Firm conclusions regarding the density dependence of 
depletions cannot be drawn from the neutral species ratios until the impact of the various other 
processes on the ionization balance (§ |5.2| ) has been evaluated. 

5.3.3. SLVvs. WLV Depletions 

We have independent estimates for the WLV and SLV depletions only for Si, Fe, Al, and 
Cr, for which the differences in depletion (WLV minus SLV) are 0.37, 0.43, 0.35, and 0.46 dex, 
respectively. While these differences are all comparable — and thus reminiscent of the correlated 
depletions found by Joseph (1988) — the small baseline in depletion precludes a definitive 
confirmation of Joseph's result. Conversely, the Si and Fe depletions are also quite consistent 
with those plotted by Sofia et al. (1994) and by Fitzpatrick (1996) for a larger range in overall 
depletion — where the difference in depletion between the two elements increases with increasing 
overall depletion. [As already noted, Savage et al. (1992) found larger differences in depletion 
toward £ Oph (component group A minus component group B) for the more severely depleted 
elements.] Toward 23 Ori, the even more severely depleted Ca appears to exhibit a somewhat 
larger difference in depletion between the WLV and SLV components (0.58 dex), but the WLV 
Ca depletion is based on considerations of the Ca ionization balance, and is therefore somewhat 
uncertain. Higher resolution observations of Ti II absorption toward 23 Ori would provide a more 
reliable estimate of the difference in depletion (WLV minus SLV) for the most severely depleted 
elements. 

5.3.4- Dust-Phase Abundances 

A number of papers have examined the "dust-phase abundances" of the more cosmically 
abundant elements in attempts to identify possible mineral constituents of the dust at different 
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levels of depletion or in different environments (e.g., Sofia et al. 1994; Fitzpatrick & Spitzer 
1997; Fitzpatrick 1997). The dust-phase abundances, defined as Aiust = A cosm i c — A gas , depend 
(obviously) on the adopted cosmic reference abundances and on the atomic data used to determine 
the gas-phase abundances. The atomic data now seem to be fairly well known for most of the 
important lines of the most abundant elements (Savage & Sembach 1996; Fitzpatrick 1997; see the 
Appendix). There are strong suggestions, however, that the total (gas + dust) abundances of at 
least some elements in the ISM may be only about 2/3 of the solar values (Sofia et al. 1992; Snow 
& Witt 1996; Meyer et al. 1997); S and (perhaps) Sn may be notable exceptions. If we use the 
solar reference abundances listed in Table |5[ then the dust-phase abundances of Mg, Si, and Fe 
are all about 30-32 per 10 6 H atoms in the SLV components, since their total abundances are very 
similar and since all three are largely depleted^. In the WLV components, Ad us t ~ 19, 25, and 32 
per 10 6 H atoms for Mg, Si, and Fe, respectively. If we assume that the total ISM abundances are 
2/3 times the solar values, then Adust ~ 17-21 per 10 6 H atoms for the three elements in the SLV 
components, and Adust ~ 6, 13, and 20 per 10 6 H atoms, respectively, in the WLV components. In 
either case, Adust (Fe+Mg) ~ 2 Adust (Si) for both WLV and SLV components, consistent with most 
of the three elements being present in the dust as olivine [(MgjFe^SiOJ. In contrast, Fitzpatrick 
(1997) found Ad us t (Fe+Mg) significantly greater than 2 Ad us t(Si) for clouds exhibiting the mildest 
depletions, and suggested that some of the Mg and Fe might be in non-silicate grains. In the SLV 
components, Ad us t(Fe) may be slightly greater than Ad us t(Mg), but Ad us t(Fe) is significantly larger 
in the WLV components. 

5.3.5. Depletions in High- Velocity Gas 

Determination of column densities for the significant ionization states of various elements in 
the HV components affords an opportunity to gauge how much dust has survived in that ionized, 
presumably shocked, high-velocity gas. As discussed above, the HV column densities of C II, N II, 
and Si II + Si III are consistent with solar relative abundances of C, N, and Si and an H II column 
density of about 5 x 10 17 cm~ 2 . The limit on N(S II) + N(S III) suggests that the total N(H II) 
is not more than a factor 1.5 higher than that value, which would imply that N(C III) < iV(C II) 
and iV(N III) < iV(N II) [assuming negligible N(S IV) and no depletion of S]. Furthermore, C thus 
cannot be depleted by more than 0.2 dex [compared to the 0.4 dex found for the predominantly 
neutral gas in sightlines characterized by a wide range in /(H2) (Sofia et al. 1997)]. Likewise, Si 
cannot be depleted by more than 0.2 dex. On the other hand, Al does appear to be depleted, 
by 0.4 to 0.6 dex — so that some dust has apparently survived. If Si is not depleted, then the 
depleted Al might be in the form of oxides (i.e., instead of in silicates) in the remaining dust. The 
uncertainties in the Si abundance would allow the dust-phase Al to be incorporated in silicate 
grains, however. More extensive data for £ Ori (including column densities for C III, N III, S III, 
and Fe III) suggest that Si may be depleted by about 0.3 dex, while Al and Fe may be depleted by 
about 0.8 dex, in the high- velocity gas along that line of sight (Welty et al. 1999c). 



We note that the relative gas-phase abundances of Mg in the WLV and SLV components were assumed equal to 
those for Si, in view of the similar average depletions of the two elements for the total line of sight and of the similar 
relative SLV-WLV abundances found for Al, Cr, and Fe. 
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5.4. High Ions 



The relatively low resolution obtained for spectra of the high ions C IV, N V, Si III, Si IV, and 
S III makes the interpretation of the low-velocity absorption from those species rather uncertain 
(§§ [O] and The observed N(C IV)/iV(N V) = 2.6±0.2 is slightly lower than the median value 
4.0±2.4 (rms deviation) found by Sembach, Savage, & Tripp (1997) for 23 lines of sight with data 
from WE and/or GHRS. The observed N(C TV)/N(Si IV) = 8.3±0.7, however, is significantly 
higher than the corresponding median value 3.8±1.9 for 31 lines of sight. Three lines of sight in the 
Sembach et al. compilation have comparably high N(C IV)/A r (Si IV) ratios, but only one (a Vir) 
also has a low N(C IV)/iV(N V) ratio. Higher resolution spectra obtained with the GHRS echelle 
have revealed complex structure in the profiles of these high ions in several lines of sight — for 
example, relatively narrow components (at similar velocities) in Si IV and C IV, intermingled with 
broader components seen in N V and C IV (e.g., Savage, Sembach, & Cardelli 1994; Fitzpatrick & 
Spitzer 1997). Some of the narrower high ion components may be associated with components 
seen at similar velocities in lower ionization stages. While the absorption from C IV, N V, Si IV, 
and S III toward 23 Ori appears to be symmetric and featureless (in the G160M spectra), slight 
differences in central velocity among the four species (Table 12) suggest that multiple components, 
arising from different locations and due to different physical processes, may be present. We briefly 
consider several possible origins for these high ions. 



• Hot gas in collisional ionization equilibrium (CIE) at log(T) ~ 6.2 K (as is thought to 
fill the Orion-Eridanus bubble) would include a small percentage of C IV and N V (e.g., 
Sutherland & Dopita 1993). The observed 6-values (23-24 km s _1 ), however, indicate that 
those ions arise in gas with log(T) <^ 5.7 K. In CIE, the peak abundances of C IV and N V 
occur at log(T) ~ 5.0 K and 5.25 K, respectively (Sutherland & Dopita 1993). The peak 
abundances in non-equilibrium situations, however, can occur at either higher T (ionization 
lagging behind heating) or lower T (recombination slower than cooling). The calculations of 
Schmutzler & Tscharnuter (1993), for example, indicate maximum abundances of C IV and 
N V in isochorically (i.e., at constant density) cooling gas at log(T) ~ 4.25 K. If the Al III 
component at —6 km s _1 is associated with the higher ions, its width implies that log(T) <^ 
4.8 K. 

• In principle, some of the observed Al III, S III, Si III, and Si IV could be produced by 
photoionization. The observed N(C IV)/V r (Si IV) ratio is larger than any of the values 
predicted by Gruenwald & Viegas (1992) for spherical H II regions, but could be consistent 
with the 04 star models if Si were slightly depleted. It is unlikely, however, that any N V 
could be due to photoionization, due to the large (77.5 eV) ionization potential of N IV. 

• Borkowski, Balbus, & Fristrom (1990) described calculations of the high ion populations 
produced in conduction/evaporation fronts at the interface between hot (10 6 K) and cool 
(<J 10 4 K) gas — as might be expected, for example, between the hot, diffuse gas within 
the Orion-Eridanus bubble and either the walls of the bubble or any cool, neutral clouds 
remaining within the bubble. The observed N(C IV)/iV(N V) ratio is consistent with 
the values predicted for times 10 5 yr t <J 10 6 yr, if the magnetic field is essentially 
perpendicular to the interface, or for brief periods slightly before 10 5 yr for fields at smaller 
angles relative to the interface. The predicted ratios are too large for earlier times (for 
any field orientation) and too small for later times. The observed N(C IV)/iV(Si IV) ratio 
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is generally lower than the predicted values, but adding photoionization of Si III to the 
models may remove that discrepancy. For t ~ 10 5 -10 6 yr, the models predict somewhat 
lower overall column densities than those observed, as well as somewhat narrower absorption 
line profiles [dominated by thermal broadening at log(T) ~ 5.0-5.5 K]. Consistent column 
densities and line widths might be achieved if several such interface regions, at somewhat 
different velocities, were present along the line of sight. 

• Slavin & Cox (1992) calculated the long-term evolution of the structure and high ion 
content of spherical supernova remnant (SNR) bubbles expanding into a uniform, low 
density medium, with and without the effects of magnetic pressure. The inclusion of a 5 /xG 
magnetic field significantly affects the structure of the bubble — producing a thicker, less 
compressed shell and a smaller interior volume of hot gas. The predicted average column 
densities for C IV, N V, and Si IV exhibit roughly constant ratios for times t ~ 0.4-5.0 x 10 6 
yr. Within that time range, the predicted iV(C IV)/iV(N V) ratio is typically ~ 2 (slightly 
smaller than the observed ratio), and N(C IV) /N (Si IV) is ~ 15 (a factor ~ 2 larger than 
observed) . Photoionization of Si III by radiation from within the remnant would increase the 
predicted iV(Si IV), however. The models predict comparable widths for the C IV and Si IV 
lines, but the Si IV lines toward 23 Ori appear to be somewhat narrower. The predicted 

C IV line widths are also smaller than those observed. 

• Slavin, Shull, Sz Begelman (1993) modelled the steady-state behavior of "intermediate 
temperature" gas [log(T) ~ 5.0-5.5 K] produced in turbulent mixing layers formed at the 
boundaries between hot and cool gas. The predicted N(C IV)/iV(N V) ratios are typically 
greater than 10 for any transverse velocity 25 km s _1 <J v t <^ 100 km s _1 between the hot 
and cool gas — much higher than the value observed toward 23 Ori. 

While none of the above models can completely reproduce the high-ion column density ratios 
and line widths observed toward 23 Ori, the observed high ions may plausibly be attributed to a 
combination of a SNR bubble (Orion-Eridanus) and one or more conductive interfaces (between 
the hot bubble gas and cooler gas either within the bubble and/or in the shell — perhaps the 
WLV and/or the IV gas, which are at similar velocities to the high ions). Some additional 
photoionization of Si III, reasonable in either case, would be needed to produce the observed 
Si IV. Higher resolution spectra of the high ions would allow more definitive tests. It would also 
be interesting to explore models in which multiple supernovae re-heat and further expand existing 
cavities (as is likely to have happened in the Orion region), and to predict the abundances of 
additional observed species such as Al III, Si III, and S III (as well as others, such as C III, N III, 
and S IV, which will be observable with FUSE). 



6. Conclusions / Summary 

We have used spectra obtained with the Hubble Space Telescope GHRS together with 
high-resolution optical spectra and UV spectra from Copernicus to study the abundances and 
physical conditions in the diffuse interstellar clouds seen along the line of sight to the star 23 Ori. 
High-resolution (0.3-1.5 km s _1 ) optical spectra of Ca I, Ca II, Na I, and K I were used to define 
the underlying component structure for the neutral gas that is unresolved in the UV spectra, 
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enabling more accurate abundances and physical conditions to be derived for individual clouds 
from the numerous transitions found in the UV. 

Multiple absorption components are present for each of several distinct types of gas, which 
are characterized by different relative abundance and depletion patterns and physical conditions. 

• Strong low-velocity (SLV) absorption, due to cool, moderately dense neutral gas and 
representing about 92% of the total iV(H I), is seen for various neutral and singly ionized 
species at +20 km s _1 <J v © ^ +26 km s" 1 . For these neutral SLV components, most 
typically severely depleted species are less depleted by factors of 2-4, compared to the "cold, 
dense cloud" pattern found, for example, in the main components toward ( Oph. For the 
blend of the two strongest SLV components, the relative populations of the J=0 and J=l 
rotational levels of H2 imply T ~ 100 K; the C I fine structure excitation then implies a 
thermal pressure log(nnT) ~ 3.1 cm _3 K, and thus wh ~ 10 — 15 cm" 3 and a total thickness 
of 12-16 pc. The average SLV electron density, derived independently from nine pairs of 
neutral and singly ionized species, ranges from about 0.04 cm" 3 to about 0.95 cm" 3 ; we 
adopt n e = 0.15+0.05 cm -3 . The relatively large n e /nn ~ 0.01 implies some ionization of 
hydrogen in these predominantly neutral components. 

• Weaker low-velocity (WLV) absorption, probably due to warmer neutral gas, is seen 
primarily for various singly ionized species at km s _1 <J v ® <J +30 km s" 1 . The depletions 
in the neutral WLV gas are typically less severe by another factor of 2-3 than in the SLV 
gas, and are somewhat similar to the "warm cloud" pattern seen in lines of sight with low 
reddening, low mean density, and/or low molecular fraction. If T ~ 3000 K for the WLV 
components, then we have log(nHT) ~ 4.7-4.8 cm~ 3 K, «h ~ 15-20 cm -3 , n e ~ 0.2 cm" 3 , 
n e /nu ~ 0.01, and a total thickness of 0.7-0.9 pc. 

• Weak intermediate-velocity (IV) absorption from ionized gas is seen in C II, Mg II, and Si III 
at —43 km s" 1 <; v® — 4 km s" 1 . This gas is denser and slightly more ionized than the 
high velocity gas, though the depletions may be comparable, and appears to arise in thin 
filaments or sheets (thicknesses <; 0.001 pc to 0.04 pc). 

• Absorption from a number of singly and doubly ionized species, perhaps due to a radiative 
shock ("Orion's Cloak", Cowie et al. 1979), is seen at -108 km s" 1 <; v® <; -83 km s" 1 . 
While the depletions in the ionized high-velocity (HV) components are uncertain, due to 
unobserved ionization stages, Al is probably depleted by a factor ~ 3, even at cloud velocities 
in excess of 100 km s" 1 ; C and Si are depleted by no more than 0.2 dex. The individual 
ionized high-velocity components typically have T ~ 8000+2000 K, n e = wh ~ 0.4-0.5 cm" 3 , 
thermal pressure log(2n e T) ~ 3.7-4.0 cm" 3 K, and thicknesses of order 0.1 pc. 

• No obvious absorption is discerned from a circumstellar H II region around this Bl V star, 
though none would be expected if 23 Ori is located within the hot, ionized Orion-Eridanus 
bubble. 

• Weak absorption from the more highly ionized species S III, C IV, Si IV, and N V is seen 
in G160M spectra, centered at —5 km s" 1 <; v® <^ +6 km s" 1 . The relative abundances of 
the high ions are consistent with an origin in a supernova remnant bubble (Orion-Eridanus), 
with additional contributions from interface regions between hot and cool gas within the 
bubble. 
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The large range in n e derived from different pairs of neutral and singly ionized species in 
the SLV gas suggests that additional processes besides simple photoionization and radiative 
recombination affect the ionization balance. Charge exchange with protons may reduce the 
abundances of S I, Mn I, and Fe I (though the relevant rates generally are not well known at T ~ 
100 K). Dissociative recombination of CH + may help to enhance C I. A slightly harder far-UV 
field in the Orion region could selectively reduce (somewhat) the abundances of neutral species 
with ionization potentials near that of hydrogen. There is no apparent relationship between n e 
and the elemental depletion determined from the dominant singly ionized species. The apparent 
systematic enhancement of Ca I found for many lines of sight remains unexplained. Previous 
determinations of n e for other lines of sight, typically based on limited data and limited knowledge 
of the potentially relevant physical processes, must be viewed as very uncertain. 

We have briefly considered several possible explanations for the large apparent fractional 
ionization (n e /nu ~ 0.01) in the SLV gas. Such a fractional ionization would seem to require 
an ionization rate Ch two orders of magnitude higher than the value typically adopted for the 
local ISM. An enhanced flux of X-rays, which appears to be present in the Orion region, could 
in principle ionize some of the hydrogen in the SLV clouds, but more quantitative estimates of 
the potential ionization rate are needed. Alternatively, if the SLV components are located at the 
boundary of the Orion-Eridanus bubble, mixing with adjacent ionized gas (estimated to have 
of order one-tenth the column density) might yield the apparent high n e . A third possibility is 
that charge exchange between singly ionized atomic species and large molecules (if present at a 
fractional abundance of order 10~ 7 -10~ 6 ) could enhance the abundances of the neutral species 
- lowering n e and Ch by factors of up to 10 and 100, respectively. A smaller n e would be more 
consistent with the values derived from the C II fine-structure excitation equilibrium. In any 
case, n e /nn is higher than the value 2 x 10 -4 that would obtain if the electrons were due only to 
photoionization of carbon. 

Comparisons of the SLV depletions and kh with those found for the strong "component B" 
blend toward C Oph (where the densities are higher by at least a factor 10 and the depletions 
of many elements are more severe by factors of 2-4) hint at a possible relationship between 
depletion and local density for relatively cold clouds. Similar data for individual clouds along 
other lines of sight are needed to confirm any such relationship, however, as the SLV depletions 
may instead reflect the history and location of the SLV components in a swept-up shell. While 
ratios of some trace neutral species might seem to suggest more severe depletions for some 
elements in the denser parts of the SLV clouds, uncertainties as to the importance of various 
processes which may affect the ionization balance of individual elements preclude firm conclusions. 
Considerations of the ionization balance in the WLV and ( Oph A components suggests that Ca 
may generally be significantly more depleted than Fe in warm, low density gas. In both the SLV 
and WLV clouds, the dust-phase abundance of Mg+Fe is roughly twice that of Si, whether solar 
or B-star abundances are adopted as the reference. The dust-phase abundance of Fe appears 
to be comparable to that of Mg in the SLV components, but is significantly larger in the WLV 
components. The apparent residual depletion of Al in the high-velocity gas (and similar data for 
the HV gas in several other lines of sight) provides constraints on the efficiency of dust destruction 
in shocks. 



Higher spectral resolution data for some of the species observed with G160M and/or 
Copernicus would enable some of the remaining uncertainties regarding individual cloud properties 
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to be addressed: 

• New spectra of H2 would yield more accurate column densities and temperatures for the 
individual SLV components. 

• Spectra of Mg I A2853, C I A1656, N I, O I, Si II* A1265, and Ti II would yield a better 
characterization of the individual WLV components (ionization state, temperature, density, 
electron fraction, and depletion); the C I data might also provide constraints on the relative 
bh in the two strongest SLV components. 

• Spectra of C III, N II, N III, Si III, and Fe III would yield more accurate abundances and 
depletions for individual ionized components at low, intermediate, and high velocities. 

• Higher resolution spectra of the higher ions C IV, N V, Si III, Si IV, and S III would yield 
a clearer picture of the relationships among those ions and of the processes responsible for 
their creation. 

Observations of Ha emission (e.g., with WHAM) should provide useful constraints on the ionized 
gas at low velocities. Finally, determinations of charge exchange rates at low temperatures are 
needed to gauge the possible effects of charge exchange processes on the ionization balance — to 
better constrain the ionization rate required in the SLV clouds and to evaluate the possibility of 
enhanced depletions of some elements at higher densities. 

Most of the observations reported here were obtained with the NASA/ESA Hubble Space 
Telescope at the Space Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. We thank D. 
Doss (McDonald), J. Spyromilio (AAO), and D. Willmarth (KPNO) for assistance in setting up 
the various ground-based spectrographs, J. Salzer for providing Ti II data, and E. Fitzpatrick for 
obtaining the KPNO Na I spectrum. Support for this work was provided by NASA through grants 
GO-2251.01-87A, GO-5882.01-94A, and GO-05896.01-94A (administered by STScI) and through 
LTSA grant NAGW-4445. 

A. Atomic Data 

In general, wavelengths and oscillator strengths were taken from Morton (1991). In a number 
of cases, however, different values were adopted — in light of more recent work or in order to 
achieve a consistent set of component parameters for all lines of a given species. 

C I: For all transitions of C, N, and O (except the intersystem transitions of C II, N I, and 
O I described below), we have adopted /-values from the recent compilation by Wiese, Fuhr and 
Deters (1996). In general, the new values differ little from those in Morton (1991), but they 
represent a new standard. For C I, there are some significant differences shortward of 1200 A, 
where Morton depended on observed interstellar lines or calculations by Kurucz & Peytremann 
(1975). We have had to retain those older calculations for the lines at 1138.6, 1140.4, and 1158.7 
A. Our data for 23 Ori suggest that the /-values for C I A1276.5, C I* A1157.8, A1189.2, and 
A1193.7 may be too high, by about 0.4, 0.5, 0.4, and 0.2 dex, respectively. The /-values for the 
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A1276.5, A1193.7, and A1157.8 lines therefore appear to be more consistent with the values derived 
by Zsargo, Federman, & Cardelli (1997) from interstellar curves of growth. 

C II: Fang et al. (1993) measured the lifetime of the upper level of the important intersystem 
transition at 2325.405 A and used a branching fraction from theory to find / = 5.78 x 1CT 8 . 

N I: The /-values for the intersystem lines at 1159.817 and 1160.937 A are from the 
calculations of Hibbert et al. (1985). 

O I: The mean of the lifetimes by Johnson (1972), Wells & Zipf (1974), Nowak, Borst, & 
Fricke (1978), and Mason (1990) for the intersystem transition at 1355.598 A is 180±7 ^sec. 
Biemont & Zeippen (1992) calculated 200 /xsec, and a branching fraction of 0.755 for the A1356 
transition. We prefer to use this branching fraction with the experimental mean to derive / = 
1.16 x 10~ 6 , compared with 1.25 x 10~ 6 in Morton (1991). 

Mg II: Fleming et al. (1998) have published a careful calculation for the weak doublet at 1240 
A which addresses the severe cancellation of terms that had troubled previous attempts. For the 
multiplet, they find / = 8.3 x 10~ 4 , which in LS coupling gives / (1239.9) = 5.54 x 10~ 4 and 
/(1240.4) = 2.77 x 10~ 4 . We have adopted these values, pending experimental confirmation. 
Fitzpatrick (1997) compared these lines with the damped profiles at 2796 and 2802 A toward ( 
Oph, £ Per, and HD 93521 by using a detailed model for the complex component structure for each 
line of sight, and concluded that the values for the weaker doublet should be /(1239) = 6.4 x 10~~ 4 
and /(1240) = 3.2 x 10~ 4 , close to the adopted values. Our similar comparisons for 23 Ori and 1 
Sco are consistent with Fitzpatrick's results. Fitzpatrick (1997) has discussed the implications of 
the revised Mg II /-values for studies of interstellar electron densities and dust composition. 

Si I: O'Brian & Lawler (1991) combined their accurate lifetimes with branching fractions 
from Smith et al. (1987) to obtain / = 0.210, with a l-a error of 0.022 dex, for the transition 
at 2515.0725 A. Another of their lifetimes implies that / (1845) = 0.207, but it suffers from 
considerable uncertainty due to unknown infrared decays. The /-value for the A1562 line appears 
to be too high by at least 0.4 dex. 

Si II: Dufton et al. (1992) performed a careful calculation which resulted in /(1808) = 0.0020, 
with an estimated uncertainty of 25%. Then Bergeson & Lawler (1993) measured a lifetime for the 
upper level which gives / = 0.00208, with a l-a error of 10%. We have adopted the latter value. 
Spitzer & Fitzpatrick (1993) revised the /-values for the lines at 1304 and 1526 A by comparing 
their profiles towards HD 93521 with that of A1808 and their adopted /(1808) = 0.0020. The 
increase to 0.00208 results in /(1304) = 0.0894 and /(1526) = 0.114. For the latter transition, 
however, we prefer /(1526) = 0.127 from the measurements of Schectman, Povolny, & Curtis 
(1998). 

S I: Biemont et al. (1998) and Beideck et al. (1994) have listed new /-values for a number of 
lines based on measured lifetimes and either experimental or theoretical branching fractions. The 
1296.2 line appears to be stronger than expected, and may be contaminated by a residual detector 
feature. 

CI I: Theoretical calculations by Biemont, Gebarowski & Zeippen (1994) support the 
measurements of Schectman et al. (1993), which we have adopted, so that /(1347) = 0.153. The 
/-value for the 1167.1 line appears to be too high by at least 0.3 dex. 
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Cr II: Bergeson k Lawler (1993) have derived /(2056) = 0.106 and /(2062) = 0.0778 from 
measured lifetimes and branching fractions. 

Mn II: We derived new values /(1162) = 0.024, /(1163) = 0.018, and /(1164) = 0.011 by 
comparison with the multiplet near 1197 A towards 23 Ori and the assumption of LS coupling 
ratios for A1164, which is blended with Ge II. The /-values for the latter multiplet depend on 
a curve-of-growth fitting of these interstellar lines with the stronger multiplet near 2600 A by 
Lugger et al. (1982). Our values are close to those which Allen, Snow, & Jenkins (1990) derived 
in a similar way for a Sco, when allowance is made for differences in the /-values adopted for the 
A1197 multiplet. 

Fe II: Improved /-values based on measured lifetimes and branching fractions are now 
available for many transitions, including /(2250) = 0.00182 from Bergeson, Mullman, & Lawler 
(1994) and /(2344) = 0.114, /(2374) = 0.0313 from Bergeson et al. (1996). Mullman, Sakai, & 
Lawler (1997) have extended these results to /(1608) = 0.058 by absorption ratios. The values 
/(1142) = 0.00247 and /(2368) = 6.28 x 10 -5 are from a comparison of interstellar lines by 
Cardelli &: Savage (1995), who tied their scale to the laboratory measurements for A2250 and 
A2260. The factor 2 decrease in /(1142) leaves uncertain the quoted /-values for the lines at 
1143.2 and 1144.9 A in the same multiplet. Cardelli & Savage also derived /(1611) = 0.00102, but 
our comparison with other lines toward 23 Ori and 1 Sco implies /(1611) ~ 0.0013, roughly 30% 
larger. This difference for /(1611) suggests that /(2368) should be higher as well. 

Co II: Mullman, Cooper, & Lawler (1998) have combined experimental lifetimes and branching 
fractions to obtain /(2059) = 0.00751 and /(2026) = 0.0125. Mullman et al. (1998) have extended 
these data with absorption measurements to obtain /(1553) = 0.0116. We still must depend on 
the calculations of Kurucz for the lines at 1187.4 and 1547.9 A. 

Ni II: Zsargo &: Federman (1998) adopted /(1412) = 0.00665 from the calculations of Kurucz 
and used observations of interstellar absorption lines to estimate /(1370) = 0.144, /(1393) = 
0.0189, /(1710) = 0.0666, and /(1742) = 0.0776 (and others that we have not observed toward 23 
Ori). Recent accurate measurements of lifetimes and branching fractions by Fedchak & Lawler 
(1999) give /(1710) = 0.0348 and /(1742) = 0.0419, indicating that the values of Kurucz and 
of Zsargo & Federman may be too large by 0.275 dex. Pending measurements of the shorter 
wavelength lines by Fedchak & Lawler, we have decreased all the values from Kurucz or Zsargo & 
Federman in Table |3| by 0.275 dex. The new /-values imply that the depletion of Ni is generally 
very similar to that of Fe (Table |5|) — so that the [Ni/Fe] ratio is not useful for determining 
depletions in QSO absorption-line systems (confirming the conjecture of Lu et al. 1996; see also 
Welty et al. 1997). 

Zn II: The rest wavelengths of the A2026 and A2062 lines seem to be inconsistent by ~ 1 km 
s _1 for a number of lines of sight; both lines were observed together with lines of other species 
(Mg I and Cr II), so the inconsistency cannot be attributed solely to uncertainties in the GHRS 
velocity calibration. 

Ge II: These /-values are based on new calculations by Biemont, Morton, & Quinet (1998), 
which supercede the values listed by Morton (1991). In the 1991 paper /(1237) appeared 
inadvertently too small by a factor 10. 

Kr I: These /-values are from recent accurate measurements by Lang et al. (1998). 
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Sn II: Migdalek (1976) derived /(1400) = 0.798 from a relativistic calculation including 
core polarization. If his relative /-values are used to scale the beam-foil lifetime of Andersen & 
Lindgard (1977), then /(1400) = 1.03, as listed in Table 3. A new experiment is needed before a 
definitive /-value can be assigned. 

Further data on wavelengths and oscillator strengths for the resonance lines of elements Ge 
and heavier are available in a new compilation by Morton (1999). 
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Fig. 1. — High-resolution optical spectra of 23 Ori, showing the strong and weak low velocity gas (0 
km s _1 <J v <J 30 km s" 1 ). The tick marks denote the components derived from fits to the various 
optical and UV spectra; some components are numbered above the KPNO Ca II spectrum. Note 
the expanded vertical scale for Ca I and the molecular species. The vertical dotted lines separate 
the IV, WLV, and SLV components. The triangle at +16.4 km s -1 denotes the zero point for LSR 
velocities; the * at +37 km s _1 denotes the stellar radial velocity. 
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Fig. 2.— Selected UV spectra of 23 Ori (most obtained with the GHRS ECH-A or ECH-B at 
resolutions of about 3.5 km s _1 ), showing primarily the strong and weak low-velocity gas. The tick 
marks denote the components derived from the profile fits; some components are numbered above 
the C II and Fe II A2374 spectra. Note the expanded vertical scale for Si I and Fe I. Absorption due 
to other species is noted above several spectra. The vertical dotted lines separate the IV, WLV, 
and SLV components. 
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Fig. 3. — Selected UV spectra of 23 Ori (most obtained with the GHRS echelle at resolutions of 
about 3.5 km s -1 or with the G160M grating at resolutions of 13-20 km s^ 1 ), covering a wider 
velocity range — to show the high- and intermediate- velocity components and the absorption from 
more highly ionized species. The tick marks denote the components derived from the profile fits; 
some components are numbered above the N II and Si II spectra. Note the expanded vertical scale 
for some weak lines. Absorption due to other species is noted above several spectra. The vertical 
dotted lines separate the HV, IV, WLV, and SLV components. The dotted line superposed on the 
S III spectrum (lower right) represents the absorption due to C I (SLV) and Si II (HV, WLV, SLV). 
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Fig. 4. — Schematic diagram of the various components seen toward 23 Ori. The asterisks denote 
discernible components (with corresponding component numbers); the solid lines show the range 
of absorption for most lines; the dotted lines show the range of absorption due to the wings of the 
strongest lines. The weak (WLV) and strong (SLV) low velocity gas is primarily neutral, while the 
gas at v <; km s _1 (IV, HV) is primarily ionized. 
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Fig. 5. — Relative abundances for individual WLV and SLV components, as derived from the profile 
fits. Uncertainties are la; limits are 2a; some components are numbered in the Zn II/S II panel at 
lower left. The horizontal dotted lines give the average values for the WLV and SLV components. 
At the right, the solar system (meteoritic) ratio (0) and the values found for Galactic cold (C), 
warm (W), and halo (H) clouds are noted. 
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Fig. 6. — Depletions for the weak (WLV) and strong (SLV) low velocity gas (left, for each element), 
ordered by elemental condensation temperature (Wasson 1985). The Galactic warm and cold cloud 
depletions (dotted and solid lines) have been adapted and updated from Jenkins (1987) (see Table 
5); the cold cloud depletion of sulfur is assumed here to be —0.2 dex. The £ Oph component (blend) 
A and B depletions (right; see § 3.2 for definitions) are from references listed in the text, adjusted 
for our adopted /-values and solar abundances. The SLV gas toward 23 Ori seems less depleted (for 
most moderately to heavily depleted elements), by a factor 2-4, compared to the typical Galactic 
cold cloud or ( Oph comp B pattern. The WLV depletions are broadly consistent with the Galactic 
warm cloud and £ Oph comp A pattern. 
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Fig. 7. — Empirical curves of growth for trace neutral species (squares), dominant, little-depleted 
species (triangles), and dominant, depleted species (circles) in the H I gas toward 23 Ori. The 
equivalent widths and oscillator strengths are given in Table 3; the column densities are given in 
Table 7; only detected, unblended lines are used. Solid lines are predicted curves based on adopted 
component structures for Na I, Zn II, and Fe II. The strongest lines of Mg II and C II exhibit 
both damping wings and intermediate- velocity absorption, and thus lie slightly above the predicted 
curves. Uncertainties (Tcr) are not plotted where they are less than ±0.1 dex. Several lines have 
apparently discrepant /-values, as discussed in the Appendix. 



-56- 



Table 1. 23 Ori — Optical Spectra 



Date 


Facility 


FWHM 


S/N 


Species Observed 


Ref 






(km s- 1 ) 








1978 Jan 


McD 2.7m + PEPSIOS 


2.0 


45 


Na I 


1 


1992 Sept 


KPNO coude feed + cam 5 


2.5-3.0 


90-180 


K I, Ca II, Ca I, CN, CH, CH+ 




1993 Oct 


KPNO coude feed + cam 6 


1.2-1.4 


65-110 


K I, Ca II, Ca I, CN, CH+ 


2 


1994 Feb 


AAT + UHRF 


0.32 


30 


Ca II 


2 


1995 Oct 


KPNO coude feed + cam 6 


1.3-1.4 


60-125 


Ca II, Ca I, CN, CH, CH+ 


2 


1995 Oct 


McD 2.7m coude 


0.56 


60 


K I 


3 


1996 Apr 


KPNO coude feed + cam 6 


1.5 


95 


Na I 


4 



References. — 1 = Hobbs (1978b); 2 = Wclty ct al. (1996); 3 = Welty et al. (1999d) ; 4 = Fitzpatrick & 
Welty (1999) 
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Table 2. 23 Ori — GHRS Spectra 



Order 


Sct a 


^min Amax 


FWHM b 


S/N 


Species c 






(A) 


(km s- 1 ) 








Echellc Spectra 


A-45 


III 


1257.5-1264.2 


3.35 


65-70 


S II, Si II, Fc II, C I, CI*, CI** 


A-42 


III 


1330.4-1337.8 


3.50 


70-75 


C II, C II* 


B-32 


I 


1738.8-1750.3 


3.61 


30-35 


Ni II, Mg I 


B-31 


I 


1801.9-1813.0 


3.54 


35 


S I, Si II, Ni II 


B-31 


I 


1821.9-1833.6 


3.39 


40 


Mg I 


B-30 


II 


1852.7-1865.1 


3.62 


55-75 


Al III 


B-28 


I 


2019.8-2031.5 


3.36 


35 


Zn II, Co II, Mg I, Cr II 


B-27 


I 


2052.6-2065.7 


3.67 


40 


Cr II, Zn II 


B-25 


II 


2242.8-2257.2 


3.49 


65 


Fe II 


B-24 


I 


2317.9-2332.4 


3.62 


125 


C II], Ni I 


B-24 


I 


2334.4-2346.8 


3.33 


35 


Fe II (spybal) 


B-24 


I 


2364.0-2377.5 


3.33 


40-45 


Fe II, Al I 


B-23 


II 


2476.9-2491.7 


3.27 


85 


Fe I 


B-22 


II 


2512.0-2527.8 


3.71 


100-110 


Si I, Fc I, Co I 


B-20 


II 


2792.0-2808.0 


3.55 


85 


Mg II, Mn I 


G160M Spectra 




II 


1114.3-1152.8 


20.0 


10-65 


C I, Fc II, P II, N I, CO 




II 


1156.8-1194.1 


19.2 


100-145 


C I, Si II, Si II*, N I, Kr I, Ge II, S III, CI I, Mn II, Co II 




II 


1229.3-1267.6 


17.8 


100-110 


C I, N V, Mg II, Kr I, Ge II, Si II, Si II*, S I, S II, Fe II 




11 


1270.2-1308.4 


17.2 


110-125 


C I, O I, Si II, P II, S I 




II 


1325.7-1363.7 


16.3 


90-110 


C I, CI I, Cu II, C II, C II*, O I, CO 




11 


1378.0-1415.8 


15.6 


70-110 


Sn II, S IV, Ni II, P I, S I, Ga II, CO 




II 


1503.2-1538.8 


14.1 


30 


P II, Si II, Si II* (spybal) 




II 


1534.6-1571.9 


13.8 


100-110 


C IV, C I, Si I, Fc II, Co II 




II 


1588.5-1625.7 


13.3 


90-110 


Fe II, Ge II, Zn I 




II 


1645.1-1681.9 


12.8 


100-125 


C I, Mg I, Al II, P I 



a Set I (program 3993) was observed on 1992 April 14-15; set II (program 2251) was observed on 1992 September 23; 
set III (program 5896) was observed on 1996 January 13. 

b For echelle spectra, resolution = 78000 + 6.4(mA — 55000), where m is order number. For G160M spectra, resolution 
ranges from A/0.076A at short wavelengths to A/0.071A at long wavelengths. 

c Potentially observable; not all were detected (see Tables 3 and 4). 



Table 3. 23 Ori - 
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1560.7090 


P 


1 


1656.2672 


P 


1 


1657.3792 


P 


1 


1657.9068 


P 


1 


1139.8120 




1 


1158.0186 




1 


1158.1315 







1158.4919 




-0 



580 d 


14.8±5.2 




—39.7 


136 


8.9±4.7 






206 


36.8±4.9 


B 


—36.6 


573 d 


12.5±4.1 


B 


-34.7 


197 d 


4.2±1.9 






451 d 


49.3±4.0 


B 




810 d 


18.0±1.8 




-41.2 


104 d 


39.9±1.8 


B 




254 d 


blend with 


Si II 




004 d 


5.4±1.3 




-33.7 


731 d 


blend with 


Si II 




184 d 


27.7±2.2 


B 


—35.3 


680 d 


37.6±1.4 




f+2 3) 


680 d 


36.3±1.6 




—35.4 


758 e 


14.6±1.4 




—36.3 


077 d 


75.0±1.9 


B 


-36.8 


467 d 


34.4±1.5 




—36.2 


923 d 


45.2±1.9 




-33.4 


050 d 


60.2±2.2 




—25.0 


362 


97.8±2.6 


B 


—24.5 


082 


<4.9 






603 


blend with 


C I 




390 


<4.8 






326 e 


<1.9 






849 d 


blend with C I 




337 


2.9±1.7 






736 d 


5.5±1.7 


B 


-34.9 


664 d 


blend with C I* 




529 e 


<1.0 






226 d 


2.4±2.1 






541 d 


blend with Si II 




090 d 


blend with Si II 




085 e 


8.0±2.1 


B 


—35.3 


323 d 


blend with C I* f 




581 d 


11.4±3.1 f 


B 




231 d 


10.7±0.8 




+2.3 


122 d 


9.2±0.8 




+2.2 


268 d 


11.5±0.8 




+2.2 


508 d 


1.9±1.3 






955 d 


blend with C I 




455 d 


30.4±1.9 


B 


-36.0 


026 d 


<1.3 






208 


12.9±1.3 




-36.0 


737 d 


5.8±1.7 


B 


-36.3 


936 d 


5.6±1.5 




-36.3 


452 d 


blend with C I* 




539 d 


43.4+2.0 


B 


-33.9 


328 d 


blend with C I* 




925 d 


54.6±2.4 


B 


-24.9 


449 d 


blend with C I* 




989 


41.3±2.2 




-24.8 


771 


33.7±2.5 




-24.3 


894 


35.3+2.1 




-24.1 



129 blend with C I 

375 d blend with C I 

627 d blend with C I 
549 d <1.5 
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Table 3 — Continued 



Ion 


\ a 

Avac 




log(fA) a 




Ad c 


W A (HV) 




(A) 






(mA) 


(km s- 1 ) 


(mA) 


C I** 


1189.4469 




0.602 d 


<1.2 








1189.6307 


S 


1.056 d 


1.5±1.2 








1193.2401 




1.639 d 


blend with Si II 








1193.3932 




1.018 d 


blend with Si II 








1193.6489 




-0.387 d 


blend with C I* 








1194.0635 




0.844 d 


blend with C I 








1261.4257 


HS 


1.105 d 


1.7±0.6 


+2.4 






1261.5519 


HS 


1.582 d 


5.4±0.9 


+2.4 






1277.5496 




2.005 d 


blend with C I* 








1277.7229 


S 


1.297 d 


2.1±1.3 








1277.9538 




0.019 d 


<1.3 








1279.2286 




0.684 d 


<1.3 








1280.3328 




1.260 d 


blend with C I* 








1280.8470 




0.799 d 


<1.1 








1329.5775 


P 


1.801 d 


10.8±1.9 B 


-33.7 






1329.6005 


P 


1.325 d 


blend with C I 








1561.3402 


P 


1.227 d 


blend with C I** 








1561.3667 


P 


0.048 d 


blend with C I** 








1561.4384 


P 


1.974 d 


nil \ £~\ a T~» 

24.1±2.4 B 


—24.9 


<1.8 




1657.0082 




2.236 


blend with C I 








1658.1212 


P 


1.771 


9.8±1.8 


-24.0 




C II 


1334.5323 


HP 


2.229 d 


267.5±3.7 


+2.1 


121.7+1.5 




1334.5323 




2.229 d 


251.5±3.5 




122.6+2.5 




2325.403 


H 


-3.872 d 


<0.6 






C II* 


1335.6627 


H 


1.229 


11 1 ■ j 1 /~1 11^ 

blend with C II* 








1335.7077 


HP 


2.183 d 


175.3±2.7 B 


+1.7 


9.3+2.3 




1335.7077 




2.183 d 


170.0±2.6 B 




11.2+2.8 


C IV 


1548.195 


P 


2.470 


37.6±4.1 




<2.6 




1550.770 


I' 


2.169 


20.7±3.3 




<2.6 


N I 


1134.1653 




1.237 d 


116.3±17.9 


-39.5 






1134.4149 




1.526 d 


109.1±22.6 


-38.1 






1134.9803 




1.691 d 


128.4±21.8 


-37.3 






1159.817 


P 


—2.006 


3.4±1.8 








1160.937 


S 


-2.553 


1.7±1.2 








1199.5496 




2.193 d 


161.3±2.8 f B 




<2.5 




1200.2233 




2.015 d 


138.3±2.3 f 


— 13.8 






1200.7098 




1.713 d 


132.3±2.5 f 


-13.9 




N II 


1083.990 


P 


2.096 d 


I57.5±4.5 s 




41.5+4.1 


N V 


1238.821 


P 


2.286 


7.l±l.l 








1242.804 


P 


1.986 


5.3±l.6 






U 1 


1QOO 1 COK 

1302. IboO 




1.831 


1 or n 1 A 


—33.8 


<1.7 




lore C^n r 7 r 7 
1355. 59 / / 


1 


— 2.803 


3.7=bl.7 


—29.1 




Na I 


5891.5833 


HP 


3.570 


215.3±3.7 


+0.1 






5897.5581 


HP 


3.273 


183.2±4.0 


+0.1 




Mgl 


1651.1636 




0.036 


<2.1 








1658.3115 




0.163 


<1.2 








1668.4288 




0.302 


3.1±1.7 


-25.7 






1747.7937 


HP 


1.213 


12.5±1.8 


-1.0 






1827.9351 


HP 


1.651 


25. Oil. 6 


-1.1 






2026.4768 


HP 


2.369 


45.1±2.3 


-1.0 


<3.8 
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Table 3 — Continued 



Ion 






log(fA) a 




Av c 


W A (HV) 




(A) 






(mA) 


(km s- 1 ) 


(mA) 


Mg II 


1239.9253 


P 


-0.163 d 


24.6±1.3 


-36.2 






1240.3947 


P 


— 0.464 d 


17.8+1.2 


—36.5 






2796.352 


HP 


3.234 


405.8+4.2 


-8.0 


77.7+2.9 




2803.531 


HP 


2.933 


369.9±3.9 


-7.0 


44.7+3.7 


Al I 


2367.7762 


HP 


2.398 


<1.8 




<3.3 


Al II 


1670.7874 


P 


3.486 


148.2+2.6 


-22.7 


20.1+2.6 


Al III 


1854.7164 


HP 


3.017 


11.6+4.3 


(-4.5) 


2.3+1.8 




1862.7895 


HP 


2.716 


6.0+2.8 


(-4.5) 


1.4+1.2 


Si I 


1562.002 




2.769° 


<1.9 








1845.5202 


S 


2.58: d 


7.7±4.7 f 








2515.0725 


HP 


2.723 d 


6.7±1.2 


-8.7 


<1.7 


Si II 


1190.4158 




2.474 


170.0+4.1 B 


—33.9 


8.6+2.2 




1193.2897 




2.775 


236.6+5.3 B 


—33.9 


18.9+1.2 




1260.4221 


HS 


3.104 


235.2+3.2 B 


(+2.9) 


44.6+1.7 




1260.4221 




3.104 


234.7+2.6 B 


—32.2 


44.6+2.0 




1304.3702 




2.067 d 


140.3+2.2 


-32.6 


3.4+1.9 




1526.7066 




2.288 d 


175.7+10.3 








1808.0126 


HP 


0.575 d 


67.4+4.1 


-0.9 




Si II* 


1194.5002 




2.872 


blend with C I* f 








1197.3938 




2.173 


<2.0 f 








1264.7377 


P 


3.058 


4.2±1.5 








1265.0020 




2.104 


<1.9 








1533.4312 




2.306 d 


<3.8 f 






Si III 


1206.500 f 


P 


3.304 


181.1+2.7 




92.8+2.1 


Si IV 


1393.755 


P 


2.855 


10.1±2.6 




<2.4 




1402.770 


P 


2.554 


4.4+2.4 




<2.4 


P I 


1381.4760 


P 


2.038 


<1.8 








1674.5953 




1.440 


<1.7 








1679.6969 




1.616 


<1.7 






P II 


1152.818 




2.435 


46. 1+10. 8 g 








1301.874 


P 


1.351 


18.2±1.8 


—32.5 






1532.533 


P 


1.067 


12.0±3.8 f 


2.2 




S I 


1247.1602 




0.742 d 


<1.9 








1270.7804 




1.039 d 


<2.1 f 








1295.6531 


P 


2.052 d 


8.2+2.0 


—33.0 






1296.1739 




1.455 d 


6.9±2.0 


-31.9 






1303.4300 




0.769 d 


<1.4 








1401.5142 


S 


1.267 d 


2.8±1.6 








1807.3113 


HP 


2.239 d 


13.5±1.9 


-0.9 


<3.4 


S II 


1250.584 




0.834 


76.3±2.5 


-34.5 






1253.811 




1.135 


95.8±2.2 


-34.4 






1259.519 


H 


1.311 


111.3±1.3 


(+1.9) 


1.4+1.3 




1259.519 




1.311 


109.0±2.5 


-34.5 


2.1+1.5 


S III 


1190.208 




1.421 


blend with Si II 




<1.4 
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Table 3 — Continued 





\ a 






vv a 


Av c 






(A) 






(mA) 


(km s- 1 ) 


(mA) 


Cl I 


1167.1480 




1.048 e 


<1.3 








1188.7742 




1.937 


blend with C I 








1335.7258 


H 


1.521 


blend with C II 








1347.2396 


P 


2.314 d 


22.8±2.0 


-30.3 




K I 


7701.093 


HP 


3.417 


43.8±3.1 


+0.4 




Ca I 


4227.918 


HP 


3.870 


4.1±0.6 


+0.4 




Ca II 


3934.777 


HP 


3.397 


74.5+1.5 


+0.0 




Ti II 


3384.740 


P 


3.061 


5.5±0.7 h 






Cr II 


2026.269 


H 


0.980 


<2.0 








2056.254 


HP 


2.338 d 


19.2±3.6 


—0.9 






2062.234 


HP 


2.205 d 


11.9±2.8 


— 1.2 




Mn I 


2795.641 


H 


3.198 


<1.0 








2799.094 


H 


3.069 


<1.0 








2801.907 


H 


2.909 


<1.0 








average 


HP 


3.07 


0.5±0.3 






Mn II 


1162.015 




1.45: d 


5.5±1.8 


-40.1 






1163.326 




1.32: d 


3.1±1.6 


-43.1 






1164.208 




l.ll: d 


blend with Ge II 








1197.184 


P 


2.273 


21.6+2.3' 


— 14.0 






1199.391 




2.104 


11 1 " j 1 TV T Tf 

blend with N I 1 








1201.118 


P 


2.024 


12.9±2.0 f 


— 12.8 




Fe I 


2484.0211 


HP 


3.141 


5.5+1.0 


—6.0 


<2.2 




2523.6084 


HP 


2.848 


3.5±0.9 


-7.9 




Fc II 


1142.3656 


S 


0.450 


6.0+2.4 


—37.0 






1143.2260 




1.182 c 


28.7±4.0 


-34.9 






1144.9379 




2.080 e 


59.7+3.7 


-33.4 






1260.533 


H 


1.499 


blend with Si II 


+2.9 






1565.5255 




-1.569 


<1.9 








1608.4511 


P 


1.970 d 


82.8+2.4 


-25.1 


<1.9 




1611.2005 




0.321 d 


9.8±2.1 


-25.8 






2249.8768 


HP 


0.612 d 


14.7±2.2 


—5.4 






2344.214 


HP 


2.427 d 


178.3±4.6 


—2.7 






2367.5905 


H 


-0.828° 


<1.8 








2374.4612 


HP 


1.871 d 


106.5±4.0 


-1.4 


<3.5 


Co I 


2522.122 


HP 


2.762 


<0.7 






Co II 


1187.417 




1.169 


<1.1 








1547.946 




1.580 


blend with C IV 








1552.760 




1.256 d 


<1.8 








2026.407 


HP 


1.404 d 


<2.3 








2059.471 


HS 


1.189 d 


<1.6 






Ni I 


2320.744 


HP 


3.201 


<0.6 






Ni II 


1317.217 


P 


2.009 d 


7.8±3.0 f 


-4.6 






1370.132 


P 


2.020 d 


10.9±2.6 f 


-5.2 






1393.324 




1.146 d 


<2.1 








1412.866 




0.698 d 


<1.6 








1741.549 


HP 


1.863 d 


9.8±2.6 


0.0 
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Table 3 — Continued 



Ion 


\ a 

A vac 

(A) 




log(fA) a 


(mA) 


(km s- 1 ) 


W A (HV) 
(mA) 


Ni II 


1804.473 


H 


0.836 d 


<1.9 






Cu II 


1358.7730 


P 


2.713 


4.0±1.6 


-27.7 




Zn I 


1589.5610 


P 


2.287 


<2.1 






Zn II 


2026.136 
2062.664 


HP 
HP 


3.018 d 
2.717 d 


73.7±4.1 
62.0±4.5 


-0.9 
-2.2 




Ga II 


1414.402 


P 


3.406 


4.6±2.0 


-25.1 




Gc II 


1164.272 

1237.0591 

1602.4863 


P 


2.513 d 
3.183 d 
2.362 d 


3.5±2.0 B 
10.6±1.9 
2.5±2.3 


-38.0 




Kr I 


1164.8672 
1235.8380 
average 


P 


2.344 d 
2.402 d 
2.37 


<1.6 
<1.9 
<1.3 






Sn II 


1400.440 


P 


3.158 


1.6±1.4 







a Vacuum wavelengths and /-values from Morton (1991), unless noted 
otherwise. H = high-resolution spectra (FWHM < 3.7 km s~ x ). P = primary 
source for column density; S = secondary source for column density. 

b B = major contributor to blend. High-velocity absorption is not included 
(see last column). Uncertainties and limits arc 2-cr. 

c Observed velocity offset relative to Ca II. 

d Revised value which differs from Morton (1991) by 0.01 dex or more; see 
Appendix for details. 

°Possibly discrepant oscillator strength; see Appendix for details. 

f Data from program 2584, obtained from HST archive. 

g Data from Copernicus; see Cowie et al. (1979) and Jenkins et al. (1986). 

h Data from Magnani & Salzer (1991). 
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Table 4. 23 Ori — Molecular Column Densities 



Molecule 


A 

(A) 


Transition 


log(fA) a 


(mA) 


N c 
(cm- 2 ) 


b c 

(km s- 1 ) 


(km s- 1 ) 


Ref 


CH 


4300.3132 


A-X (0,0) 


1.341 


4.0±0.7 


4.9±0.4el2 


2.4±0.2 


23.1±0.1 


W99 


CH+ 


4232.548 


A-X (0,0) 


1.367 


9.4±0.9 


11.4±0.4el2 


2.7±0.1 


23.6±0.1 


W99 


CN 


3874.602 


B-X (0,0) 


2.117 


<1.1 


<2.5ell 






W99 



CO 



H 2 



1150.483 
1344.186 
1392.525 
1544.448 

1108.128 
1108.634 
1096.439 
1079.226 
1115.896 
1096.725 
1081.710 
1053.976 
1085.144 
1070.898 
1125.539 
1089.512 
1113.868 
1097.843 



B-X (0,0) 
A-X (7,0) 
A-X (6,0) 
A-X (0,0) 

0-0, R(0) 

0- 0, R(l) 

1- 0, P(2) 

2- 0, R(2) 

0- 0, P(3) 

1- 0, R(3) 

2- 0, R(3) 
4-0, R(3) 

2- 0, R(4) 

3- 0, R(4) 
0-0, P(5) 

2- 0, R(5) 

3- 1, P(l) 

4- 1, R(l) 
total 



0.716 
0.745 
1.305 
1.387 

0.283 
0.113 
0.398 
0.902 

-0.105 
0.589 
0.887 
1.175 
0.881 
1.059 

-0.087 
0.832 
1.196 
1.389 



<4 
<3 
<3 
<4 

77±15 
73±14 
32±15 
51±13 
46±15 
66±15 
75±21 
71±10 
<25 
33± 7 
21±15 
<28 
<41 
<11 



<2.3el3 
<4.4el3 
<l.lel3 
<1.0el3 

4.0±1.0cl7 
1.2t jel8 
1.0±0.5el7 

3.0±1.0el7 



4.0±2.0el5 
1.5±1.0cl5 
<4.7cl3 



2-0ig;fel8 



W99 
W99 
W99 
W99 

FJ,W99 
FJ,W99 
B,W99 
FJ 

FJ,W99 
B 
FJ 
B 

FJ,W99 
B 

FJ,W99 

FJ 

FJ 

FJ 

W99 



References. — W99 = this study; B = Bohlin et al. (1983); FJ = Frisch & Jura (1980) 

a Oscillator strengths from Morton & Noreau (1994) for CO; from Morton & Dinerstein (1976) or Frisch & Jura (1980) for 



b Uncertainties and limits are 1-a. 
Uncertainties are l-cr; limits are 2-<r. 
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Table 5. Galactic Abundance and Depletion Patterns 





Solar 





Depletions 




Element 


Abundance a 


Cold b 


Warm b 


Halo c 


C 


8.55 


-0.4 


-0.4 


(-0.4) 


N 


7.97 


-0.1 


-0.1 


(-0.1) 


O 


8.87 


-0.4 


-0.4 


(-0.4) 


Na 


6.31 


-0.6 


-0.4 


(-0.2) 


Mg 


7.58 


-1.1 


-0.5 


-0.3 


Al 


6.48 


-2.4 


-1.1 


(-0.6) 


Si 


7.55 


-1.3 


-0.4 


-0.3 


P 


5.57 


-0.5 


-0.2 


(-0.1) 


g 


7.27 




0.0 


0.0 


CI 


5.27 


-0.6 


-0.2 




K 


5.13 


-0.6 


-0.4 


(-0.2) 


Ca 


6.34 


-3.6 


-2.0 


(-0.8) 


Ti 


4.93 


-2.8 


-1.2 


-0.7 


Cr 


5.68 


-2.1 


-1.2 


-0.6 


Mn 


5.53 


-1.5 


-1.0 


-0.7 


Fe 


7.51 


-2.2 


-1.4 


-0.6 


Co 


4.91 


-2.1 


-0.9 




Ni 


6.25 


-2.2 


-1.4 


-0.6 


Cu 


4.27 


-1.4 


-1.1 




Zn 


4.65 


-0.4 


-0.2 


-0.1 


Ga 


3.13 


-1.1 


-0.6 




Ge 


3.63 


-0.8 


-0.5 




Kr 


3.23 


-0.2 


-0.2 


(-0.2) 


Sn 


2.14 


-0.1 


0.0 


( 0.0) 



a Solar system mctcoritic (CI carbonaceous 
chondritic) or photospheric (for C, N, and O) values 
from Anders & Grevesse (1989) and Grevesse & Noels 
(1993); logarithmic, with H = 12.00. 

b Updated from Jenkins (1987); see Lauroesch et al. 

(1996) and Welty ct al. (1997, 1999a). Differences 
from previous versions are due to use of meteoritic 
reference abundances (P, S, CI, Ti, Mn), new /-values 
(C, Mg, Co, Ni, Ge; see Appendix), and recent studies 
of C II [Sofia et al. (1997)], N I [Meyer et al. (1997)], 
O I [Meyer ct al. (1998)], and Kr I [Cardclli & Meyer 

(1997) ]. Warm cloud value for Ca estimated in this 
paper (§ 3.2.3). Depletions for heavy elements Cu, 
Ga, Ge, and Sn are based on values in Hobbs et al. 
(1993), Cardelli (1994), and Lambert et al. (1998.) 

c Adaptcd from Savage & Scmbach (1996); 
Fitzpatrick (1996b); and references therein. Values 
in () arc estimated. 
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Table 6. 23 Ori — H I Region (WLV, SLV) Component Structure 1 



Comp 


VQ 

(km s- 1 ) 


Na I 

5891 
5897 


K I 

7701 


Ca I 
4227 


Ca II 

3934 


Si II 

1808 


S II b 
1259 


Fe II 
2250 
2344 
2374 


Zn II 

2026 
2062 


10W 


(0.4) 


0.4±0.3cl0 
(2.75) 


<0.2cl0 


<0.3cl0 


0.6±0.2cl0 
(2.5) 


1.0±0.9el3 
(2.6) 


3.0±0.5el3 
(2.55) 


2.2±0.2el2 
(2.45) 


<0.3cll 
(2.4) 


11W 


4.9±0.2 


8.8±0.4el0 
(2.9) 


<0.2el0 


<0.3cl0 


4.9±0.2cl0 
2.7±0.3 


6.0±1.3cl3 
(2.8) 


16.0±1.0cl3 
(2.75) 


8.2±0.3cl2 
(2.65) 


3.0±0.4ell 
(2.6) 


12W 


9.9±0.8 


6.2±0.5el0 
(2.75) 


<0.2el0 


<0.3cl0 


2.2±0.2cl0 
(2.5) 


6.6±1.8cl3 
(2.6) 


16.0±1.0cl3 
(2.55) 


5.0±0.4el2 
(2.45) 


2.7±0.5ell 
(2.4) 


13W 


13.0±0.3 


1.0±0.4el0 

(2.5) 


<0.2el0 


<0.3cl0 


3.8±0.2cl0 
2.0±0.3 


6.0±2.6el3 
(2.15) 


7.0±0.6cl3 
(2.05) 


3.1±0.5el2 
(1.9) 


1.3±0.5ell 
(1.85) 


14W 


(16.05) 


1.4±0.4el0 
(2.5) 


<0.2cl0 


<0.3cl0 


1.4±0.2cl0 
(2.0) 


5.3±4.1el3 
(2.15) 


4.5±0.5cl3 
(2.05) 


5.9±0.9cl2 
(1.9) 


1.3±0.5ell 
(1.85) 


15W 


18.4±0.6 


2.0±1.0cl0 
(1.85) 


<0.2cl0 


<0.3cl0 


2.3±0.3cl0 
(1.5) 


(7.3cl3) 
(1.65) 


9.0±5.0cl3 
(1.6) 


11.6±3.7el2 
(1.35) 


1.5±1.0cll 
(1.35) 


16S 


(20.1) 


7.5±2.0cll 
1.3±0.3 


1.2±0.5el0 
(1.3) 


<0.3cl0 


5.5±0.5el0 
(1.5) 


33.±14.el3 
(1.5) 


(45.el3) 
(1.5) 


20.±8.el2 
(1.5) 


8.0±3.4cll 
(1.5) 


17S 


21.65±0.1 


(100.±40.cll) 
0.9±0.15 


13.7±0.4cl0 
0.9±0.1 


1.0±0.1cl0 
(0.9) 


39.±5.el0 
1.05±0.1 


51.±17.cl3 
(1.05) 


(201.cl3) 
(1.05) 


54.±4.3cl2 
(1.05) 


37.±8.9cll 
(1.05) 


18S 


23.80±0.1 


(115.±50.cll) 
1.35±0.2 


16.1±0.4el0 
1.35±0.15 


0.6±0.2el0 
(1.35) 


63.±6.el0 
1.50±0.2 


133.±43.cl3 
(1.5) 


(630.cl3) 
(1.5) 


109.±8.cl2 
(1.5) 


114.±28.cll 
(1.5) 


19S 


26.2±0.3 


3.1±0.6cll 
0.8±0.3 


0.4±0.2el0 
(0.8) 


<0.3cl0 


2.8±0.3cl0 
(1.0) 


17.2±7.2cl3 
(1.0) 


(73.cl3) 
(1.0) 


16.1±4.6el2 
(1.0) 


13.3±5.3ell 
(1.0) 


20W 


28.05±0.3 


l.lil.OelO 
(1.7) 


<0.2cl0 


<0.3cl0 


1.8±0.2cl0 
(1.3) 


11.7±2.2el3 
(1.5) 


18.0±7.0cl3 
(1.4) 


6.3±1.6el2 
(1.15) 


4.1±1.7cll 
(1.1) 


a Entries 


are N (cm" 


~ 2 ) and b (km s~ 


r ). Values in 


parentheses were fixed in the fitting. Uncertainties are 1- 


-<j; limits are 2- 


<t. Components 



(16-19) are SLV; components (11-15;20) are WLV. N 17 /N 18 fixed in fitting UV lines of Si II, S II, Fe II, and Zn II (§ 3.1.2). 

b Total column density, and thus SLV column densities, determined by using curve of growth (§ 3.1.3). Uncertainties for WLV components 
reflect variations for changing velocity offset by ±0.3 km s and SLV column densities by ±0.3 dex. 
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Table 7. 23 Ori — HI Region (SLV, WLV) Column Densities and Depletions 21 



Ion 



log[N(tot)] 
(cm -2 ) 



Dtot 



log [N (WLV)] 
(cm -2 ) 



log[N(SLV)] 

(cm -2 ) 



D SLV Method b 



H I 

H 2 



2D 74+ - 08 
18.301°- 



(19.61) 



(20.71) 



18.30 



+0.11 
-0.07 



3,5 
3,5 



c 


I 


14.93 


c 


I* 


13.99 


c 


J** 


13.18 


c 


Itot 


14.99 


c 


II 


16.98 


c 


II* 


14.6 



+0.05 
-0.05 
+0.08 
-0.11 
+0.05 
-0.07 
+0.05 
-0.05 
+0.06 
-0.08 



-0.31 



12 56+° 11 
12 52+°' 06 
12 11+ 009 
12 90+ 07 
(15.85) 
14.3-14.6 



(-0.31) 



14 OQ+0.05 

„ 7 +0.08 
-0.09 
io i 5 +0.05 

14 98 +0 ' 05 

lt.3O_ Q5 

(16.95) 
14.2-15.0 



(-0.31) 



2,4,5 

2,4,5 

2,4,5 

2,4,5 

3 

2 



N I 
O I 

Na I 



16.64loTo 
17.32_°;}_ 



13.36±0;i 2 



-0.15 
-0.35 



(15.51) 
(16.19) 



ii OO+0.04 
1±.A4_ Q4 



(-0.15) 
(-0.35) 



(16.61) 
(17.29) 



13.35t°;lJ 



(-0.15) 
(-0.35) 



Mgl 
Mg II 

Al I 
Al II 



13.81_°;°_ 
15.68.^ 

<11.53 

io 55+O.17 



-0.64 



-1.67 



11 95+0.20 

(14.88) (-0.31) 



12.73 



+0.02 
-0.02 



-1.36 



13.80l°;°2 ■ ■ ' 2 .5 
(15.61) (-0.68) 3,5 



<11.53 



13.48 



+0.17 



-1.71 



Si I 
Si II 
Si II* 



o.oi 

0.02 
0.07 
0.08 

11.6+0.1 



11.79 
15.44 



-0.85 



14.64t° ;°? 
11.5+0.1 



-0.52 



11.79 
15.37 



+0.01 
-0.02 
+0.08 
0.11 

<10.7 



-0.89 



2,5 
1,5 



P I 
P II 



<12.11 



14.04 



+0.04 
-0.04 



-0.27 



(12.91) (-0.27) 



<12.11 
(14.01) 



(-0.27) 5 



S I 
S II 

CI I 

K I 



12.81 
16.01 



+0.03 
0.04 
+0.10 



13 sn+°' 02 

11 5n+ ' 01 
ii.ou_ 01 



0.00 



14.87 



+0.05 



0.05 
(11.22) 



-0.01 



12.81 



+0.03 



(15.98) 
(13.29) 



11 50+ ' 01 
ii.ou_ 01 



2,5 

(0.00) 1,5 
5 
1 



Ca I 
Ca II 

Ti II 

Cr II 



10 20+ ' 06 

12 10+ ' 03 
iz.iu_ 03 



ii 9 o+0.03 

19 7 o+0.01 
^■'•^-O.Ol 



11.23±g;g| (-2.43) 
-2.44 (10.48) (-2.06) 



-1.69 12.00_°'o3 -1.29 



10.20_°, ; « 6 5 
12.04_°;°» 



(11.15) 



12.64+.°,;°,! 



-3.01 1 
(-2.49) 5 
-1.75 2 



Mn I 
Mn II 

Fe I 
Fe II 



10.23_£2 



11.34_°;» 
14.38±g;gi 



-1.12 (12.35) (-0.79) 



-1.87 13.63t°'o4 -1.49 



10 23+ ' 20 
iu.zo_ 40 

(13.08) 

11 34+ 03 
14 qn+0.03 



(-1.16) 5 



-1.92 1 



Co I 
Co II 



<10.73 

<12.71 <-0.95 



<10.73 
<12.71 



<-0.92 



Ni I 
Ni II 

Cu II 

Zn I 
Zn II 



<10.26 
12.98+°,;!° 



11.; 



hO.03 
-0.03 



<11.89 

, o 97 +0.06 
io.z(_ Qg 



-2.01 (12.23) (-1.63) 

-1.09 (11.06) (-0.76) 

-0.12 12.14+ ,; , 7 , (-0.12) 



<10.26 
(12.90) 

(11.79) 

<11.89 
13 24+ 07 



(-2.06) 2,5 
(-1.13) 5 

(-0.12) 1 
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Table 7 — Continued 



Ion 


log[N(tot)] 

(cm -2 ) 


D tot 


log[N(WLV)] 

(cm -2 ) 


Dwlv 


log[N(SLV)] 
(cm- 2 ) 


Dslv 


Mcthod b 


Ga II 


11 19+ 003 


-0.68 


(10.39) 


(-0.35) 


(11.12) 


(-0.72) 


5 


Ge II 


11 89+ - 04 


-0.48 


(10.76) 


(-0.48) 


(11.86) 


(-0.48) 


5 


Kr I 


<11.72 


<-0.25 






<11.72 


<-0.22 




Sn II 


10-85±g;« 


-0.03 


( 9.72) 


(-0.03) 


(10.82) 


(-0.03) 


5 



a WLV and SLV entries in parentheses are estimates based on fits to echelle spectra of "similar" 
species (§§ 3.1.3 and 3.1.4). Uncertainties are 1-cr; limits arc 2-<r. 

b l = fit individual components; 2 = scaled high-resolution (WLV, SLV); 3 = total from damping 
wings; 4 = scaled low-resolution (WLV, SLV); 5 = scaled low-resolution (total) 
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Table 8. 23 Ori — Relative Electron Densities and Abundances (H I Regions) a 



Element 


log[N(X I)] 


log[N(X II)] 


r/a h 


n c c 


[X/C]„ d 


[X/C]: c 




(cm -2 ) 


(cm -2 ) 


(cm -3 ) 


(cm -3 ) 






SLV Components 








(WJ1/D) 


(WJ1/D) 




(WJ1) 


C 


14.98 


16.95 


24./40. 


0.26/0.43 






Na 


13.35 


(14.52) 


2.2/2.6 


(0.15)/(0.18) 


(-0.19) 


-0.43 


Mg 


13.80 


15.61 


13.5/12.1 


0.21/0.19 


-0.37 


-0.46 


Al 


<11.53 


13.48 


253./299. 


<2.8/<3.4 


-1.40 


<-0.36 


Si 


11.79 


15.37 


357./368. 


0.09/0.10 


-0.58 


-1.02 


P 


<12.11 


14.01 


97./143. 


<1.2/<1.8 


+0.04 


<+0.72 


S 


12.81 


15.98 


74./114. 


0.05/0.08 


+0.31 


-0.40 


K 


11.50 


(13.34) 


10.1/11.4 


(0.15)/(0.16) 


(-0.19) 


-0.44 


Ca° 


10.20 


12.04 


66./77. 


0.95/1.11 


-2.70 


-2.13 


Mn 


10.23 


13.08 


26. /28. 


0.04/0.04 


-0.85 


-1.70 


Fe 


11.34 


14.30 


37./37. 


0.04/0.04 


-1.61 


-2.41 


Ni 


<10.26 


12.90 


41./40. 


<0.09/<0.09 


-1.75 


<-2.19 


WLV Components 








(100K/3000K) 


(100K/3000K) 




(100K/3000K) 


C 


12.90 


(15.85) 


24./198. 


(0.027)/(0.22) 




... / ... 


Na 


11.32 


(13.62) 


2.2/22. 


(0.011)/(0.11) 


(+0.01) 


-0.38/-0.29 


Mg 


11.95 


(14.88) 


13.5/137. 


(0.016)/(0.16) 


( 0.00) 


-0.23/-0.14 


Ca+ 


11.23 


13.95+D Ca 


0.05/0.45 


9.5e-5<5 -1 /8.6e-4<5 -1 


Dca-Dc 


-2.14/-2.10 



a Values in parentheses arc estimates based on assumed depletions (Na II and K II) or component structure 
(WLV C II and Mg II). For WLV Ca+, entries are for Ca II and Ca III, and D is the logarithm of the Ca depiction 
(S). 

b F/a from Aldrovandi & Pcquignot (1973; 1974), Pcquignot & Aldrovandi (1986), and Shull & van Steenberg 
(1982). For SLV, entries are for WJ1 and Draine radiation fields, both with a at T=100 K; for WLV, entries are 
for WJ1 field, for a at T=100 K and 3000 K. 

C SLV electron densities are average values over 4 SLV components [dominated by two strongest components 
(17 and 18), for which n e (17) ~ 2.6 n e (18) (§ 3.1.2)]. WLV electron densities are average values over 7 WLV 
components. 

d log(X II/C II) — (Ax — Ac), where Ax is log(mctcoritic abundance). 

e Rclativc abundance calculated from neutral species ratio, adjusted for solar abundances and (T/a) ratios (for 
WJ1 field). 



Table 9. 23 Ori - 



Physical Conditions in SLV and WLV Gas 



Property - - - SLV - - - - - - WLV - - - 



N(R I) (cm- 2 ) 


20.71 


19.61 


N(C 1*)/N(C I)tot 


0.10±0.02 


0.40±0.08 


N(C l**)/N(C I)tot 


0.015±0.003 


0.16±0.04 


T(K) 


100«o 


(100) a 


(3000) a 


log(n H T) (cm" 3 K) 


3.1±0.1 


4.3±0.2 


4.7-4.8 


71h (cm -3 ) 


10-15 


200 


15-20 


n e (cm- 3 ) 


0.15±0.05 (0.015) b 


0.02 


0.20 


n e /n H 


0.01 0.001 


0.0001 


0.01 


Nn/n H (pc) 


12-16 


0.1 


0.7-0.9 



log[7V(C II*)] obs (cm" 2 ) 14.2-15.0 14.3-14.6 

log[JV(C II*)]p r ed (cm- 2 ) 15.3 14.6 14.5 14.25 

log[iV(Si II*)] obs (cm- 2 ) <10.7 -11.5 

log[AT(Si II*)] pr ed (cm- 2 ) 10.8 9.9 10.0 11.3 



a Assumed representative temperatures for cold or warm gas. 

b Assumed lower n e — if, for example, charge exchange with large 
molecules is significant. 
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Table 10. 23 Ori — H II Region (HV, IV) Component Structure 11 



Comp 


v © 


C II 


C II* 


Mg II 


Al III 


Si II 




(km s- 1 ) 


1334 


1335 


2796 


1854 


1260 










2803 


1862 





1H 



-108.0±0.4 



0.7±0.1el3 
(3.5) 



0.6±0.2el2 
(3.5) 



0.7±0.1ell 
(2.5) 



<2.0cl0 
(2.4) 



l.l±0.2ell 

(2.35) 



2H 



-100.5±0.1 



6.4±0.4el3 
4.2±0.2 



1.8±0.2el2 
(4.2) 



11.5±0.3ell 
2.9±0.1 



9.6±2.3el0 
(2.9) 



14.6±0.5ell 
2.9±0.2 



3H 



-93.3±0.1 



3.2±0.1el3 
3.4±0.3 



1.0±0.2el2 
(3-4) 



5.9±0.3ell 
2.8±0.2 



4.9±2.3el0 
(2.8) 



11.6±0.5ell 
2.75±0.15 



4H 



51 



61 



71 



81 



-83.3±0.1 



-42.6±0.3 



-32.4±0.4 



-19.8±0.8 



-11.8±0.4 



5.6±0.3el3 
3.8±0.1 

3.2±0.2cl2 
(3.6) 

3.1±0.2cl2 
(4.5) 

3.3±0.9cl2 
(4.5) 

3.4±0.2cl3 
(6.0) 



1.4±0.2el2 

(3.8) 

<0.4el2 



0.7±0.2el2 
(4.5) 

1.0±0.2el2 
(4.5) 

(0.6el2) 



6.0±0.1ell 
2.2±0.1 

<1.3cl0 
(2.35) 

<1.4cl0 
(3.2) 

3.5±0.2el0 
(3.2) 

5.4±0.3ell 
(5.5) 



2.2±2.0el0 
(2.2) 



16.4±0.3ell 
2.25±0.07 

9.2±1.6cl0 
(2.35) 

2.8±1.7cl0 
(3.2) 

1.0±0.4ell 
(3-2) 

12.1±0.6ell 

(5.5) 



91 



-4.3±0.3 



6.7±1.3el3 
(4.5) 



(1.6el2) 
(4.5) 



.8±0.7ell 
(3.0) 



5.1±0.4cll 
6.2±0.6 



21.5±1.8ell 
(3.0) 



WLV ■■■ (7.1el5) 3.0±0.8el4 (6.8el4) ••■ (4.4el4) 

SLV ■■■ (8.9cl6) (9.0el4) (3.6cl5) 1.9±0.3ell (2.3el5) 

21W (+30.0) 0.6±0.1el4 •■• 0.2±0.1cl3 ■•■ 0.3±0.1cl3 

(3.5) ■■• (2.5) ••■ (2.5) 



a Entries are N (cm -2 ) and b (km s~ 1 ). HV, IV values in parentheses were fixed in the fitting. 
Uncertainties are 1-ct; limits arc 2-cr. Components (1—4) are HV; components (5—9) are IV. WLV 
and SLV values in parentheses refer to predominantly neutral gas (Table 7). 
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Table 11. 23 Ori — H II Region Column Densities and Depletions' 1 



Ion 


log[N(HV)] 


log[N(IV)] 


log[N(high)] 


log[N(LV)] b Method c 




(cm ) 


(cm- 2 ) 


(cm -2 ) 


(cm -2 ) 



HI <15.9 <16.2 ■■• (20.74) 

H II 17.7 to 17.9 17.7 to 18.0 >18.2 <19.3 

CI <12.04 <12.04 ■■• (14.99) 

C II 14.20±°;ol 14.05^0.^6 ' ' ' (16.98) 1 

CII* 12.67+°;g* 12.891°;^ ■■• 14.6-15.1 1 

C IV <11.75 • •• 13.00l°;°J "• 6 

D(C) -0.2to+0.0 d -0.5 to -0.2 d ■■• (-0.31) 

NI <12.25 <12.00 •■• (16.64) 

Nil 13.73±° °* 14.1±0.2 ■■• <15.3 5 

NV <11.79 ■•• 12.58jl°;°| ■■• 6 

D(N) -0.1 to +0.1 +0.1 to +0.4 •■• (-0.19) 

OI <12.36 <12.70 •■• (17.32) 

Mgl <11.96 <11.96 •■• (13.81) 

Mgll 12.38l°;°J 12.16±g;g| •■• (15.68) 1 

All <11.80 <11.80 •■• (<11.5) 

A1II 11.69t°'°5 H.7+0.3 ■■• (13.55) 5 

A1III 11.221^ 11.7ll°;g* ■■• H.281S;™ 1.2 

D(A1) -0.6 to -0.4 -0.5 to -0.2 ■■• (-1-67) 

Si I <11.28 <11.28 •■• (11.79) 

Si II 12.64115;^ 12.55+°;o2 • ■ ' (15.44) 1 

Sill* <11.15 <11.15 •■• 11.6+0.1 6 

Si III 13.11±°° 7 7 13.2+0.2 ■■• <13.5 5 

Si IV <11.38 ■•• 12.08±g"_3l ■■• 6 

D(Si) -0.2 to +0.0 -0.3 to +0.0 ■■• (-0.85) 

SI <12.03 <12.03 •■• (12.81) 

S II 12.95l°; 2 * 13.01+ ( ( ;;[* • ■ • (16.01) 2 

S III <12.70 ■•• 13.48±g:°7 •■• 6 

D(S) [0.00] [0.00] [0.00] (0.00) 

Fe I <10.86 <10.86 •■• (11.34) 

Fell <12.12 <12.12 •■• (14.38) 

a Uncertainites are l-<r; limits are 2-ct; values in square brackets are assumed. 

b Entries in parentheses arc for predominantly neutral (H I) components (WLV 
+ SLV). 

c l = fit individual components; 2 = scaled high- resolution (HV, IV); 5 = scaled 
low-resolution; 6 = fit low-resolution 

d Does not include possible contributions from (unobserved) C III. 
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Table 12. 23 Ori — Weak H II and High Ion Component Structure 



Ion 


A 

(A) 


JV a 
(cm -2 ) 


fe a 

(km s~ ) 


(km s ) 


Si II* b 


1264.7377 


3.6±0.4ell 


6±3 


12±1 


Al III 


1854.7164 


5.3±0.5ell 


6.3±0.8 


-5.7±0.5 






1.7±0.4ell 


4.7±1.7 


20.5±1.1 




1862.7895 


4.3±0.6cll 


5.6±1.1 


-4.2±0.7 






2.5±0.5ell 


3.4±1.1 


23.9±0.7 




both 


5.1±0.4ell 


6.2±0.6 


-6.1±0.4 






1.9±0.3ell 


4.4±1.1 


20.4±0.7 


S III 


1190.208 


3.0±0.4el3 


16±4 


5±2 


C IV 


1548.195 


9.8±0.4el2 


24±1 


3±1 




1550.770 


10.7±0.5el2 


22±2 


9±1 




average 


10.0±0.3cl2 


23±1 


6±3 


Si IV 


1393.755 


1.2±0.1el2 


18±2 


-6±1 




1402.770 


1.0±0.1el2 


13±3 


-3±2 




average 


1.2±0.1el2 


17±2 


-5±1 


N V 


1238.821 


3.5±0.2cl2 


23±2 


2±1 




1242.804 


5.1±0.5el2 


27±3 


-2±2 




average 


3.8±0.2el2 


24±2 


0±2 



a Uncertainties are 1-cr. 

arise in H I gas (§ 3.3.3). 
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Table 13. 23 Ori — H II Region (HV, IV) Component Properties' 1 



Comp 


jv(c n*)/iv(c ii) 


rp b 

± max 


n e = n ii c 


log(2n H T) c 


AT(H II) d 


r 






(K) 


(cm- 3 ) 


(cm" 3 K) 


(cm -2 ) 


(PC) 


1H 


0.085±0.031 


8700/8900/9150 


1.4±0.5 


4.4 


2.0el6 


0.005 


2H 


0.028±0.004 


12550/11950/13950 


0.47±0.09 


3.9 


1.8cl7 


0.12 


3H 


0.031±0.006 


8200/11150/12550 


0.52±0.12 


3.9 


9.0el6 


0.06 


4H 


0.025±0.004 


10250/6850/8400 


0.42±0.08 


3.8 


1.6cl7 


0.12 


51 


<0.125 


9200/ •■■ / ••• 


<2.1 


<4.5 


9.0el5 


>0.001 


61 


0.23±0.07 


14400/ ••• / ••■ 


3.9±1.2 


4.8 


8.7el5 


0.0007 


71 


0.30±0.10 




5.0±1.7 


4.9 


9.3cl5 


0.0006 


81 


<0.118 




<1.9 


<4.5 


9.6cl6 


>0.02 


91 


0.091±0.039 




1.5±0.7 


4.4 


1.9el7 


0.04 



a Components (1-4) are HV; components (5-9) are IV. 
b From 6-valucs for C II, Mg II, and Si II. 
c Assuming T = 8000±2000 K. 

d Prom N(C II), assuming solar photosphcric abundance, no C III, and no depletion. 



